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TECHNICAL NOTE 3030 

A METHOD FOR CALCULATING THE SUBSONIC STEADY- STATE 
LOADING ON AN AIRPLANE WITH A WING OF ARBITRARY 
PLAN FORM AND STIFFNESSl 
By W. L. Gray and K. M. Schenk 

SUMMARY 


A method for con^juting the steady-state span load distribution on 
an elastic airplane wing for specified airplane weights and load factors 
is given. The method is based on a modification of the Weissinger 
L-method and applies at subcritical Mach numbers . It includes the 
effects of external stores and fuselage on the spanwise loading. Modi- 
fications are outlined for treating tail-boom and tailless airplane con- 
figurations and for calculating the divergence dynamic pressure of a 
swept wing with a large external store. A method is also outlined for 
reducing wind-tunnel data to obtain effective aerodynamic coefficients 
which are free of model flexibility effects. The effects of Mach num- 
ber can readily be evaluated from the aerodynamic coefficients thus 
obtained. 


INTRODUCTION 


The inclusion of the effects of flexibility in the solution of the 
spanwise airload distribution applied to a wing of arbitrary plan form 
and stiffness distribution has increased the complexity of analysis over 
that for a rigid wing. The methods that are available at the present 
time are generally concerned with the calculation of loading on an iso- 
lated flexible wing rather than the more practical case not only where 
the effects of fuselage and nacelles on the spanwise loading must be 
taken into account but also where the total lift on each of the major 
components must be considered simultaneously in order to determine the 
wing loading at a specified load factor. A method for including such 
effects without recourse to iterative procedures for steady-state flight 
conditions and subcritical Mach numbers is presented in this paper. The 

%ased on Boeing Airplane Company Document No. D-1062k, "a Matrix 
Solution for the Subsonic Steady State Aeroelastic Loading on Airplanes 
by W. L. Gray and K. M. Schenk, June 1, 1951* Acknowledgement is made 
to Messrs. Paul W. Harper and John B. Garvin of the NACA for extensive 
work in editing and revising this document. 
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equations are derived so that the spanwise airload distribution can be 
expressed in matrix form in terms of influence coefficients for aero- 
dynamic induction and structural deflection in a manner similar to that 
employed in reference 1. 

The basic method is outlined in the body of the paper. Included 
in appendixes are details of the various derivations, the expansion of 
the basic equations to include fuselage interference and store load 
effects, the modifications for tail-boom and tailless configurations, 
a method for determining divergence dynamic pressures for swept wings 
with large external stores, a method for reducing wind-tunnel data to 
obtain effective aerodynamic coefficients which are free of model flexi 
bility effects, and a method for obtaining compressibility corrections. 

SYMBOLS 

The following symbols appear in the body of this report. Addi- 
tional symbols which appear only in the appendixes are defined as they 
are introduced. 


b 

wing span, in. 


fuselage lift coefficient in presence of wing at cxj. = 0, 
Lp/qS 


fuselage pitching -moment coefficient about c/4 point in 
presence of wing at ot^ = 0, 1^/qSc 


rate of change of fuselage lift coefficient with o^, 
per radian 

(‘^)a 

rate of change of fuselage pitching-moment coefficient 
with o^, per radian 

c 

local chord parallel to plane of symmetry, in. 

c 

wing mean geometric chord, in. 

X 

wing section pitching-moment coefficient 

El 

effective value of product of modulus of elasticity and 
wing section beam bending moment of inertia, lb-in. ^ 

GJ 

effective value of product of shear modiilus of elasticity 
and wing section polar moment of inertia, lb-in. 2 
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h semispan of horseshoe vortex, in. 

Ix airplane rolling moment of inertia, lb-in. ^ 

ly airplane pitching moment of inertia, lb-in. ^ 


Lp 


fuselage lift in presence of wing. 



lb 


I 


% 





^si 

”>o 


n 

P-T 

P 

P 

<1 

r 

S 


wing running lift per inch of span perpendicular to plane 
of symmetry, Ib/in. 

fuselage pitching moment about c/k point in presence of 
wing, + (^Cjiip^^a^qSc, in- lb 

airplane rolling moment^ in- lb 

elastic-airplane rolling moment caused by unit aileron 
deflection^ in- lb 

elastic-aii^lane daniping moment in roll caused by unit 
wing-tip helix angle (pb/2V)j^, in-lb 

elastic-airplane rolling moment caused by unit rolling 
angular acceleration acting on the wing distributed 
inertia, in- lb 

elastic-airplane rolling moment caused by unit spoiler 
deflection, in-lb 

two-dimensional lift-curve slope per radian, including com- 
pressibility effects, for sections parallel to plane of 
symmetry 

airplane load factor, positive when inertia loads are downward 

balancing tail load, positive upward, lb 

airplane rolling velocity, radians/sec 

airplane rolling angular acceleration, radians /sec^ 

dynamic pressure, Ib/sq in., pV^/2 

radial distance from vortex core, in. 

wing area, sq in. 
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true free-stream velocity, in. /sec 

airplane gross weight, lb 

wash velocity induced by line vortex at perpendicular dis- 
tance r from vortex line, positive for downwash, in. /sec 

downwash angle at three-quarter-chord point induced by 
vortex system representing wing and its spanwise lift 
distribution 

streamwise distance from pitch reference axis to bound 
portion of horseshoe vortex, positive when vortex is 
to rear of pitch reference axis, in. 

streamwise distance from pitch reference axis to Zfh line, 
positive when 'c/h line is to rear of pitch reference 
axis, in. 

streamwise distance from pitch reference axis to airplane 
center of gravity, positive when center of gravity is to 
rear of pitch reference axis, in. 

streamwise distance from pitch reference axis to center of 
pressure of balancing tail load, positive when center of 
pressure is to rear of pitch reference axis, in. 

lateral distance from wing center line, in. 

final angle of attack of section zero-lift line with respect 
to local free-stream direction, cxp -h ag + ag, radians 

(see fig. l) 

change in section angle of attack due to aerodynamic twists 
and due to all structural twists associated with a flexi- 
ble wing which are not accounted for by the Ug term, 
radians (see fig. l) 

angle of attack of root-section zero-lift line, radians 
(see fig. l) 

change in section angle of attack due to wing lift distri- 
bution acting on a flexible wing = 0 for a rigid wingj, 

radians (see fig. l) 

strength of line vortex, in.^/sec 


aileron deflection, radians 
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5s 

n 

e 

A 

% 


spoiler deflection 
dimensionless spanwlse station, 

airplane pitching angular acceleration, positive for nose up, 
radians/sec^ 

local sweep angle of elastic axis, radians 

equivalent local sweep angle including congress ibility 
effects, radians 


P mass density of ambient atmosphere, slugs/cu in. 

= 0 . 11 ^ 4-679 X 10“^ Ib-sec^/in.^ at standard sea- 
level conditions^ 

Matrix notation: 


[] 

LJ 


square matrix, elements of which are designated by use of 

subscripts; for example, element a^. is in ith row and 
Jth column 


row matrix 



column matrix 


diagonal matrix, which is a square matrix in which all 
elements are zeros except those on the principal 
diagonal a 22 ^ ^53^ * * * ®'nn 






aerodynamic-induction or downwash matrix in which ele- 
ments relate downwash angle at station i to 

unit running lift at station j on wing 

elasticity matrix in which elements aij relate changes 
in streamwise angle of attack at station i to unit 
running lift at station j on wing 

fuselage image-vortex matrix relating image downwash 
effects at station control points to unit rvinning lifts 
(see appendix E) 

fuselage "overvelocity" matrix (see appendix E) 


6 


NACA TN 3030 


Fil identity matrix; that is, diagonal matrix in which 

^ diagonal elements are equal to unity 

PRESENTATION OF METHOD 


In this section of the report the basic equations necessary to the 
method are outlined and discussed in a general way. Details of the 
derivations are contained in the various appendixes. 


Assumptions 

In the development of the method certain assumptions that are 
common to airfoil theory apply, namely: 

(1) The flow is potential; that is, boundary-layer effects, separa- 
tion, and compressibility shocks are absent or negligible. 

(2) The wing thickness is small. 

(3) A stagnation point exists at the wing trailing edge. 

(li-) The angles of attack a are small so that tan a sin a ^ a 
(where a is measured in radians) and cos a 1. 

(5) All drag- load effects except those due to nacelles and stores 
are neglected entirely in determining the deformations of the wing used 
in obtaining the equilibrium spanwise airload distribution. 

With regard to the structure the following assumptions are made: 

(1) Caniber changes arising from twisting and bending of the wing 
are neglected entirely. 

(2) The elastic twist of the control surface is the same as that of 
the adjoining wing structure. 

(3) The angles of structural deflection 9 are small so that 

tan 9 ** sin e ^ e (where 9 is measured in radians) and cos 9 1. 

(1<-) Although the angle -of -attack changes, including those due to 
bending and torsional deformations of the wing, are acco\xnted for in 
the determination of the equilibria spanwise airload distribution on 
the wing, this final airload distribution is applied to the geometry of 
the undeflected wing in computing the bending and torsional moments. 
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Basic Eqviations 


Symmetrical flight conditions . - The fundamental problem involved is 
the development of a series of eq^uations which relate the spanwise lift 
distribution for an arbitrary wing plan form in a given flight condition 
to the properties and attitudes of the individual sections that form the 
wing. 

If the two-dimensional wing is considered first, the following 
relationships for lift and downwash behind an airfoil are available 
from most standard textbooks on aerodynamics : 


The circulation F is taken to be such that, at a specified dis- 
tance r behind the lifting line, the resultant of the downwash veloc- 
ity Wj^ and the flight velocity V is parallel to the section zero-lift 

linej that is, no flow exists normal to the zero-lift line at this point. 
Then, 



( 1 ) 


I = m^pa^ ^ c 


( 2 ) 



( 3 ) 





and from equations (l) and ( 2 ), 



( 5 ) 


Substituting equation (5) into equation (3) results in 
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or 


VTj. ^ IDq c/2 

V~ ~ irt r 


(7) 


In order to 


satisfy equation the expression 


i“o c/2 . 

— — — in equa- 
2jt r 


tion ( 7 ) must be equal to 1.0. Since the theoretical section two- 
dimensional lift-curve slope is equal to 2%, r must equal c/2, which 
is the distance between the lifting line and the three-quarter-chord 
point . 


In the development of the method presented in this report, equa- 
tion (7) is always used in the form 



(8) 


This simplification requires that the section lift-curve slope iBq be 
the two-dimensional value (i.e., the value of the lift-curve slope for 
an unswept two-dimensional wing) and that the location of the downwash 
control point D (see fig. 2) be one-half of the local strearawise chord 
to the rear of the quarter-chord point, or at 5c A* 

The essential difference between a two-dimensional wing and a wing 
of finite aspect ratio arises from the nonuniform spanwise loading which 
produces the trailing vortices of the finite-aspect-ratio wing. The 
equations presented thus far are considered to apply to the finite- 
aspect-ratio wing when the effects of all the vortices, both bound and 
trailing, have been taken into account. 

Equation (8) in matrix form is 


o 



This matrix relation represents a series of equations, each applicable 
to a particular station on the semispan of the wing. The values of 



every one of which is affected by the bo\md and trailing 
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vortices at all of the wing stations, can be eval\iated from 



which, in combination with equation (l), resiilts in 




( 10 ) 


The 



matrix in these equations is the aerodynamic -induction or 


downwash matrix which is derived in appendix A. 


Combining equations (9) and (lO) gives 


or 



The series of equations represented by the matrix equation (l2) expresses, 
for any given dynamic pressure, the relationship between the spanwise 

variation of running lift the final section angles of attack , 

and the spanwise variation of the two-dimensional section lift-c\irve 
o 

slope The effects of wing plan-form geometry are accounted for 

thr ough the elements of the matrix. The section lift-curve slope 

is expressed in the general form mQ rather than 2« to permit substi- 
tution of actual values when available from scaled-model tests or to 
permit correction for compressibility effects as described in appendix A. 
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The final angle -of -attack variation across the span can be 

considered to be con 5 )Osed of three essential parts (see fig. l) 




H = H" N/ M 

X i/ 

For a wing free of external stores j the angle of attack caused 

by structural deflection of a flexible wing due to the section lifts 
acting at the section aerodynamic centers is linearly related to the 

matrix by an expression derived in appendix B as 


(15) 


H = PH 0} 


(It) 


The wing gp.nmetry and stiffness are accounted for in the structural- 

deflection matrix This matrix is based on loadings associated 

with stations which are psirallel to the airplane plane of symmetry. In 
a swept wing, however, the structure is usually arranged such that the 
wing boxes are formed between stations approximately perpendicular to 
the elastic axis. In order to obtain a closer representation of the 

loadings and deflections on this actml wing, an [S2^ matrix was also 
derived (see appendix B) and may be substituted for the matrix in 

equation (l^J-) when desirable. 

The matrix of eqmtlon ( 15 ) } described in detail in 

appendix C, is composed of built-in twist, apparent or aerodynamic twists 
such as those due to interference, control deflection, and ang\ilar veloc- 
ities, and all structural twists of an elastic wing which are not accounted 


for in the 




matrix. 


Although equation (l2) is perfectly general, it is not useful in 
the form given for determining the lift distribution on a flexible wing 

since a component of the matrix is itself a function of the lift. 

If therefore expressed as in equation (l4), equations (l2) 
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and ( 13 ) may be combined so as to express the load distribution on a 
flexible wing in terms of wing root angle of attack and any combination 


of the 




twists as 


0 

1 


hqm^ 




<r 

( 15 ) 


or 


0 

1 


4(111.0 


5l| - pi] 


{'} ■ H * {“g} 


(16) 


Design conditions, however, are usually specified in terms of gross 
weight and load factor rather than root angle of attack. The inclusion 
of these two additional independent variables requires two additional 
equations. By considering that the airplane must be in equilibrium as 
regards vertical forces and pitching moments, the two additional equa- 
tions may be written as .c. m . k 

It- ^e,ic,r\c>n^ * t. 4^ 


Lahlp} 


+ P, 


T 


nW = 0 


( 17 ) 


for equilibrium of vertical forces and 




-2[zh^ + 2q|2hc^ jcji^ - PrpXrp + nWx^ = 0 


‘^■3- 

( 18 ) 


for equilibrium of pitching moments about the pitch axis. 

Equations (16), ( l 7 )> and (I8) are the ba sic equations for a 
flexible-wing airplane. They may be solved simultaneously for 

and Pip as functions of any design values of speed, gross weight, 
and load factor. Equation (16) as written applies to symmetrical flight 
conditions, but by substitution of an antisymmetrical matrix 

(see appendix A) it is then applicable to unsyrametrical flight condi- 
tions which are considered in greater detail under the following sec- 
tion heading. 
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Equations (l6), (l?), and (l8) can be altered as shown in the fol- 
lowing paragraphs to include (a) the effects of forces on the fuselage, 
(b) the effects of external stores by the method derived in appendix D, 
and (c) the interference effects on the exposed-wing loading due to the 
presence of the fuselage by the method derived in appendix E. 

In connection with forces on the fuselage (item (a)), the lift and 
moment characteristics of the fuselage in the presence of the wing are 
ass\mied to be known. The lift and pitching moment of the fuselage may 
then be written with small error as 


Lp = qS 


Cl- 


■Fo 


(^Lf) 


Or 


a 


(19) 


Mp = qSc 


Vo ^ 




( 20 ) 


in whic h fus e lage up loads and nose-up mo ments are considered positive. 
This lift and moment may be appropriately included in equations (l6), 
(17), and (18) to get the following more complete set of equations 
(see fig. 2): 7 > 

For the wing load distribution. 






[z] - {oy]+ W{o} = 


( 21 ) 


for the sunmiation of vertical forces ^ 


2[_2hJ (l) + <ls(cLp)^a^ + Ft = nW - qSC^p^ 


( 22 ) 
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and for the summation of pitching moments, 



The twist terra %|o| appearing in equation (21 ) has been included 

at this point to provide for the possibility that tail loads may enter 
the wing at some point along the span, as for a tail-boom type of con- 
figuration, for instance. This contribution is otherwise con- 

sidered to be zero. A method for handling the tail-boom type of airplane 
as well as the case of the tailless airplane is described in appendix F. 

In considering the effects of the external stores (alteration (b)), 
as in the case of the fuselage, the lift and moment characteristics of 
the stores in the presence of the wing are assumed to be known so that 
the lift and moment can be given by expressions similar to equations (l9) 
and (20) for the fuselage. 

The lift -distribution equation for the store case derived in appen- 
dix D, with the Pr term included as before for generality, is 



|agj. + qS|B| (2^) 

in which the matrices |^aJ and |bJ relate an ag type of twist of 

the elastic wing to the store lift and moment. This twist is handled in 
two parts so that the part that varies with the span loading may 

be introduced on the left-hand side of the equation in a manner parallel 


to that for the 


matrix. 


ih 
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The vertical-force and pitching -moment equations ( 22 ) and ( 25 ) are 
modified to include store lift and moment in a manner similar to that 
in which the fuselage lift and moment were previously included. In the 
notation introduced in appendix D these equations for one or more stores 
become 



for simnnation of vertical forces and 


2 [ 2 hxJ + 2 qS 




I- o -1 
1 


4qmo 


[Bl] 




- qSc^Cmp^^ar + PqiXrji = 2 q [^c^J - 



for summation of pitching moments. In these equations the subscript n 
is the store index which, for any particvilar configuration, has as many 

values as there are stores on the semispan. The matrix notation [E^J 

which gives the angle of attack of the store is defined in appendix D. 

The third and final alteration (item (c)) to be considered in con- 
nection with equations (16), (17), and (18) is the change in span load 
distribution of the exposed wing due to the presence of the fuselage. 

This interference effect is derived in appendix E. The over -all fuselage- 
interference effect is considered to be conq)osed of two parts. The first 
is that due to the image vortex system within the fuselage which is 
required to satisfy the condition of zero velocity normal to and at the 
fuselage surface. This condition is satisfied by adding a correcting 
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matrix (see appendix E) to the [S]J downwash matrix. The second 

part of the effect is the increment in vertical velocities over the 
exposed wing due to the presence of the fuselage at an angle of attack. 


This effect is calculated as an interference twist of the 




type 


and is expressed as a function of fuselage angle of attack times the 

o 

”overvelo.city" matrix 


|SqJ. (See appendix E. ) The lift distribution 


(eq. (21)) when altered to include these fuselage -interference effects 
becomes 


0 

r 1 ~ 

[2"] 4. [qTI 

4 qmo 

L y L y 







(27) 


where the elements of the 



matrix give 


the increments in vertical 


velocities along the span. 


The calculation of these fuselage effects would not he required if 
appropriate data were available from wind-tunnel tests of a scaled model 
of the subject airplane. A method of determining these and other aero- 
dynamic twists as well as the applicable values of section lift-ciirve 
slope from appropriate wind-tunnel data is given in appendix G. The 
method utilizes equation (l 2 ) to obtain aerodynamic coefficients which 
are free of model wing flexibility effects and which are therefore appli- 
cable to the full-scale airplane having a wing flexibility different from 
that of the model. 


Unsymmetrical flight conditions .- In addition to the symmetrical 
flight conditions already outlined, a number of unsymmetrical flight 
conditions axe usvially investigated in structural design. Among the 
conditions which may readily be investigated by the methods of this 
report are those which arise through the use of roll-producing devices 
such as ailerons or spoilers. The load distributions on an elastic wing 
associated with roll-control deflections may be thought of as the sxirama- 
tion of distributions from the following specific loadings: 
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(1) The symmetrical or pre-roll loading which existed prior to the 
deflection of the roll-producing device. 

(2) The incremental loading resulting from the deflection of the 
roll-producing device. 

(5) The incremental loading associated with a constant rolling 
velocity with no roll-control deflection, which is usually described 
as the damping-in- roll loading. 

(14.) The incremental loading caused hy the rolling angular accelera- 
tion. This loading results from the structural twists ^Q-g^ of 
elastic wing when the rolling acceleration p acts on the wing mass 
distribution (see appendix C). Note that the resulting aerodynamic 
rolling moment will generally be in the same direction as the applied 
rolling acceleration. 

Three flight rolling conditions will be used to outline the pro- 
cedure for determining unsymmetrical loadings on the wing. These 
rolling conditions are: 

(a) Steady roll at some specified value of wing-tip helix angle pb/2V 
with no rolling acceleration 

(b) Roll initiation resulting from the instantaneous deflection of 
the roll-control device to the angle required to obtain the specified 
value of pb/2V but with no rolling velocity 

(c) Roll termination, that is, control surfaces deflected in 
opposition to the steady rolling velocity pb/2V 

Steady roll: In a steady- roll condition the span load distribution 

for the elastic wing is given by the summation of the first three loadings 
eniimerated. 


The distribution obtained for the first, or pre-roll, loading is 
described in detail in the section Symmetrical flight conditions. 
Equations (16), (I7), and (18) or their appropriate equivalents are 
used together with the symmetrical [Sj^ matrix of appendix A. 


The distribution of loading and its associated rolling moment 


•ai 


or Mjj resTOlting from a unit antisymmetrical deflection or 5s^ 

of the roll-control device are obtained from equation (I6), the antisym- 
metrical matrix jS]^ of appendix A, and the appropriate matrices 


which give the aerodynamic and structural twists resulting from a unit 

deflection of the roll-control device. These |agj matrices can be 

written in terms of control-surface deflection by the method described 
in appendix F or by means of data from wind-tunnel tests, if available. 
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Similarly, the distrihution of loading and the dainping moment 
associated with a xmit value of the wing-tip helix angle (pb/ 2 V)^ 
obtained from eq\iation (16) and the antisymmetrlcal matrix. 


are 


The 


|agj values in this case vary linearly and antisymmetrically across the 
span from (pb/2V)2^ at one tip to -(pb/2V)2^ at the other. 


These unit load distributions associated with 60, or 60-, and 

ai Si 

(pb/av)^^ must then be scaled up or down according to the amount of con- 
trol deflection 5 a or 6g required to give the desired value of pb/ 2 V. 

The deflection required is obtained from the equation of equilibrium of 
the airplane in roll as 


M- 


‘X. 


6 


1 = 


% 


■Pi 


ai 2 V (pb/ 2 V)^ 


(28a) 


or 


5 


s 




2 V (pb/ 2 V)^ 


(28b) 


where the rolling moment produced by the control deflection balances 

the rolling moment due to damping in roll in the steady specified rolling 

condition. 


After the unit load distributions have been scaled in the manner 
just described, they may be added to the pre-roll loading to obtain the 
load distribution for the specified steady-roll condition. 

In this outline the assumption of equal and opposite roll-control 
deflections is made. If, as is more generally the case^ unequal deflec- 
tions of the ailerons or spoilers are involved, the span loading must be 
determined in a slightly different way. To illustrate the procedure, a 
spoiler deflection on only one wing may be considered to be equivalent 
to a symmetrical and an antisymmetrlcal deflection with an amplitude 
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equal to half of the spoiler deflection 5 s/ 2 . The antisyrametrical 
deflection res\ilts in a gain in lift on one wing which is exactly bal- 
anced by a loss in lift on the other wing so that a rolling moment is 
produced without a change in the over-all wing lift. The symmetrical 
deflection, however, results in a change in total wing lift, and hence 
load factor, with no rolling moment. In order to compensate for this 
change, both the vertical-force and the pitching -moment balance equa- 
tions must be introduced so that the wing load distribution associated 
with the change in Oj. required to compensate for the change in tail 
load can be determined. Adding the symmetrical and unsymmetrical loading 
distributions will give the proper lift distribution for the wing with 
a single spoiler deflection 6g. 


Differentially operated ailerons might be considered in a manner 
similar to that outlined for the spoiler. A f-urther extension to include 
the combined deflection of both ailerons and spoilers can also be made 
in the same way by making use of the proper spoiler-to-aileron gearing 
ratio . 


Roll initiation: In the roll-initiation condition where no rolling 

velocity is assumed to exist, all the listed loadings occur except that 
due to damping in roll. 


Since the control deflection ( 5 a or 63) will already be known 

from the steady-roll condition, the problem is to determine the initial 
rolling angiTia-p acceleration p for instantaneous control deflection. 
The procedure involved is first to find the wing spanwise airload dis- 
tribution and its rolling moment 1*^. due to a unit rolling angular 


acceleration 


"Pi 

The values of depend on the wing mass and 


Pi- --P1 

stiffness distributions as well as wing aerodynamics. The values are 
obtained from equation (16), or its equivalent, for antisymmetric flight. 
With the value of known, the desired angular acceleration p is 

then found by solving the following equation of motion for the airplane 
in roll: 


P = 


5a 



6 


a 


1 


lx 



(29a) 
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or 


P = 


% 


Si 

Si 


T %• 

_ ^Pl 

586.4 i>i 


(29b) 


where the moments are given in inch-pounds. The value of p obtained 
from equation (29) is then used to scale the loadings previously found 
for the unit rolling acceleration to the correct value. 


Roll termination: In the roll- termination condition, the airplane 

is assumed to be rolling with a wing-tip helix euagle pb/ 2 V and the 
roll-control device is moved abruptly in a direction such as to reduce 
the rolling velocity to zero. As in the roll-initiation condition, the 
desired airplane rolling acceleration is obtained from the equation 
expressing the equilibrium of the airload and inertia- load rolling 
moments, which for aileron control is 


pb 


P = 


2 V (pb/ 2 V)i 


ix 

386.4 




■Pi 


Pi 


The airloads on the wing are those caused by the pre-roll condition 
plus the airloads from aileron deflection, danqjing in roll, and rolling 
acceleration. The inertia loads are those arising from the pre-roll 
condition plus the effects of the rolling acceleration p. 


DISCUSSION 


The method outlined in this report not only includes several previ- 
ously omitted items which are of practical interest in the design of a 
wing for aeroelastic effects but also is sufficiently extensive in scope 
that almost any type of airplane configuration may be considered. Because 
of its length, however, the method is better adapted to the determination 
of loads on a specific airplane rather than to preliminary design studies 
of several configurations. 
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Matrix formulation of the problem has particular merit for such a 
general treatment since discontinuities in angles or masses due to either 
special aerodynamic or structtiral featiores can readily be included. It 
also permits future in 5 )rovements to be made to the details of the separate 
appendixes without essentially modifying the method which has been outlined. 
In fact, throughout the paper it is possible that the engineer would modify 
the method to suit his own needs and draw from the appendixes whatever mate- 
rial would be required to investigate the problem at hand. In such a modi- 
fication he would of necessity consider the relative merits of ease in com- 
putation against the accuracy both of the method and of the data available. 

For these reasons, only a few general guides which might be considered 
for successful application of the method are given. 

For mo equal to 2n, equation (12) will give essentially the 
same results as those given by the Weissinger L-method of reference 2 
which is valid for wings of arbitrary plan form and having flat-plate, 
circular-arc, or parabolically cambered airfoil sections (refs. 3? 
and 5). The method would be expected to give the most accurate results 
when applicable values of mo are known, such as those obtainable by 
the method of appendix G for instance, since in general the fuselage, 
flaps, and external stores will affect the applicable values of mo. 

In fact, equivalent values of any of the aerodynamic parameters as 
obtained from experimental data by the method of appendix G are pre- 
ferable to purely theoretical values and may easily be Incorporated. 

The treatment of compressibility effects used in this report, 
wherein each wing section is permitted to have its own compressibility 
correction, differs from the Prandtl-Glauert method in that the wing 
p lan form is not distorted; instead, the angles of attack are altered 
as indicated by equation (ll). The treatment adopted has the merit of 
considerable saving in time for equal or better accuracy since only 
one [§il matrix is required for all Mach nmbers. The methods of 
obtaining compressible values of mo are described in appendixes A 
and G. 

With regard to the number and selection of the horseshoe vortices 
to be used to represent span loading, it is suggested that the horse- 
shoes be chosen narrower over that portion of the span where large gradi- 
ents in loading are expected, that is, near the ends of control stirfaces, 
near large changes in sweep, and at the wing tips. At least two vortices 
should be used with each control surface and a minimum of seven per semi- 
span is suggested for a clean wing. 

With regard to the structural parameters El and GJ required, it 
may be stated that equivalent values which include the stiffnesses con- 
tributed by the leading- and trailing-edge structure should be used in 
preference to the usually conservative values enqployed in the structviral 
analysis of the wing for shear and bending stresses. 
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In the design of a simple wing-fuselage combination without external 
stores or nacelles, it is necessary to study and apply the results of 


matrices as well as the fuselage -interference effects. Effects of store 
and nacelle loads and moments on the wing load distribution are covered 
in appendix D. Appendix F outlines the modifications required to adapt 
the method to the determination of wing loads on flexible tailless 
tail-boom airplane configurations. Appendix G indicates a procedure for 
obtaining equivalent values of section lift-curve slopes, effective com- 
pressible section sweep angles, and interference twists from wind-tixnnel 
tests of models which may not he scaled correctly for flexibility. Appen- 
dix H deals with the determination of the divergence dynamic pressure of a 
swept-wing airplane with a large external store. The problem of divergence 
norms.lly does not occur with a swept wing except that the attachment of 
a large external store may cause it to diverge. The determination of 
the divergence dynamic pressure is the only case in this paper which 
requires iterative procedures. 


Boeing Airplane Company, 

Seattle, Wash., July 8 , 


only appendixes A, B, C, and E which describe the 



1955. 
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APPENDIX A 


AERODYHAMIC FUNDAMENTALS 


The 



Matrix 


The pixrpose of this appendix is to explain, in simplified form, the 
aerodyneimics involved in and the steps necessary for the computation of 

the downwash matrix fi] and to develop the correction for congress i- 
hility used in this report. 

The lift or circulation distribution vhich varies along the span 
of a ving can be visualized as resulting from a system of horseshoe 
vortices, each of which is of constant strength. Such a system of 
horseshoes is illustrated in the following sketches, in which double 
arrows are used to indicate that the sense of circulation around each 
line-vortex segment is given by the right-hand rule: 




Actual airload 
curve 

Approximation to the 
actual loading as 
given by horseshoe 
vortices 
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It is obvious from the sketches that the shape of the actual load 
distribution may be approximated to any practical degree of accuracy by 
a suitable change in the number of horseshoes, each of constant strength. 
The point of importance is that the net strength of the trailing vortex 
at any point on the span of the wing is niomerically eqml to the rate 
of change of strength of the bound vortex in the spanwise direction. 

The strength of the trailing vortices would therefore be greater for 
those portions of the wing span over which the more rapid changes in 
the spanwise airload distribution occur. 

Results of theoretical investigations have shown that little loss 
in accuracy with respect to the spanwise airload distribution will be 
entailed if: 

(1) The total strength of the chordwise system of bound vortices 
is concentrated in one bound vortex located at the local streamwise 
quarter- chord point. 

(2) The downwash angle at each vortex station across the span of 
the wing, at a point one-half of the local streamwise chord downstream 
of the bound vortex, is equal to the geometric angle of attack. Herein- 
after, this point is referred to as the downwash control point D. 

r>— a.f-1 O'f' 

The downwash angle at any such control point D is therefore the 
total induced downwash velocity at that point, normal to the plane of 
the wing and caused by the complete system of bound and trailing vortices , 
divided by the flight velocity of the wing. 

It should be mentioned that the condition described in paragraph (2) 
is true as written only for airfoils having a two-dimensional lift-c\irve 
slope equal to 2 n. As is discussed in a subsequent section of this 
appendix, the condition described in paragraph ( 2 ) is modified to the 
following form when the section two-dimensional lift-curve slope is 
different from 2 rc: 



(Al) 


The downwash velocity at a point due to a single horseshoe which 
is composed of a bound vortex and its associated pair of trailing vortices 
is known to be proportional to the strength of the circulation of that 
horseshoe and, therefore, proportional to the running lift on that por- 
tion of wing span represented by the bound vortex of that horseshoe. The 
downwash angle at any one downwash control point thus is the sum of the 
incremental downwash angles due to each of the horseshoes in the system 
of horseshoes which represent the wing and its lift distribution. 
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Within the limitations which accrue from the use of a limited nijmber 
of horseshoe vortices to represent a wing, the elements with which to 
solve the following problem are now available; 

the genmet-rv of a wi ng plan form, the angle -p f -attack varia - 
-h.-irmj and the aont.inn +Tjr.-.^iniPnqlrmHl lift mirvf> - R lQpe-^ mr 1 at 1 on across 
the span, determine the spanwise airload distribution. The unknowns 
are, of coiirse, the values of the running lift at each of several points 
on the span, A necessary condition to the determination of these unknowns 
is that as many independent eqmtions be available as there are unknowns. 
This condition can be fulfilled if the angle of attack is known at each 
of the wing stations for which the loading is to be determined. 

It is obvious that, if the strength of each bound vortex represents 
the average airload over its own portion of the wing span, good accuracy 
will be obtained if the values of the running load, as determined from 
the solution of the simultaneous equations, are considered to be valid 
at the midspan point of each bound vortex. The shape and distribution 
of the continuously varying airload curve is then obtained by fairing 
a curve through all of the points thus obtained, with the restriction 
that the loading must drop to zero at the wing tip. 

The method for determining the matrix is now ill\xstrated for 

a typical wing shown in the following sketch, which includes a system 
of horseshoe vortices and associated downwash control points: 
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In actual practice it will lisually be necessary to use a larger number 
of horseshoes to represent the wing. Four per semispan are used in 
the following presentation for convenience only. The following typical 
information is available at each spanwise station; 


The section lift-curve slopes (in two-dimensional flow) at sta- 
tions T) = r\ 2 f • • • ™o " °^1’ ' ■ ■ ’ angles of attack 

of the section zero-lift line at stations t) = 
oc*^ — ^ ^2 ^ “ * * * 


Since a linear relationship exists between the strength Fj of a 
particular horseshoe vortex j and the downwash velocity w^j at a 

particular point i on the wing plan form due to that horseshoe vortex, 
the following general eq.uation can be written: 




“ij ' 


o 1 


(A2) 


where K is a constant. A particvilar horseshoe vortex, such as 

then cavises the following values of downwash velocity at control points 1 
to 8: 


’^ll 

= KiiFi 

^51 

= Sl^’l 

1 — I 

= ^ 2 l"l 

'" 6 l 

= «^ 6 /i 

'^31 

= K^lF^ 

r— 1 

= Kyil’l 

w 

= K. r 

■w„ 

= r 

in 

41 1 

81 

81 1 


V (A3) 







Similar relationships 

exist between F2, 

. . . Fg and the 

control 

points 1 to 8, that is. 




"12 " ^12^2 

^^22 " ^22^2 

”82 ■ ''82'2 



»i3 = 

^^25 "" ^25^5 

''85 ' ’'83'’3 


(A4) 




''28 ' ’=28''8 


''88 ■ '^ 88‘'8 
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I 
I 

If the horseshoes are symmetrical with respect to the center line ' 

of the wing plan form, the following values of K are equal for the 
typical wing sketched: 1 


Kll = Kq8 

^12 “ ^87 

%5 

= ^86 

Ki 4 = Kq 5 ^ 

Kgl = Ky8 

K22 = K77 

K25 

= ^76 

K24 = 

K51 = K 58 

K32 = K87 

^55 

= K66 

K54 = 

^41 “ ^58 

K42 = 

K45 

= ^56 

K44 = K55 

>> 

K51 = K48 

K52 = K47 

^55 

= K 46 

K54 = K45 

K6i = K58 

Kg2 = K57 


= ^56 

^64 = K55 

K71 = K28 

K72 = K2Y 

K73 

= K26 

K74 = K25 

% = ^18 

^82 ^ ^17 

^83 

= ^16 

K 84 = K ^5 


Further, for symmetrical wing plan forms at zero yaw it is always 
possible to obtain a spanwise airload distribution which is unsymmetrical 
with respect to the plan-form center line as the s\im of two airload dis- 
tributions, one of which is symmetrical and the other of which is antisym- 
metrical with respect to the plan-form center line; for example. 



I 
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At each spanwise station the S component is that due to the symmetrical 
distribution of loadj the A con^jonent is that due to the antisymmetrical 
distribution, and the U component represents the algebraic sum of the 
S and A components, that is, the \ansymmetrical distribution of load 
across the span. This division offers a considerable reduction in the 
amount of work required in that, for either symmetrical or unsymmetrical 
flight conditions, airload distributions need be determined on only one- 
half of the wing, provided, of coiorse, that proper account is taken of the 
sign of the circulations existing over the other half of the wing. 

For a symmetrical distribution of airload over the span 



f2 - Fy 


(A 6 ) 


and for an antisymmetrical distribution of airload 


F 


1 


= -^8 


F 


2 ~ 




(at) 


The toted, downwash velocity at einy control point is the sum of the 
downwash-velocity contributions at that point that are induced by each 
of the horseshoe vortices in the system that represents the wing; that 
is. 


*1 - hi * “12 * ”15 * “14 + "15 "16 * "17 * "18 
h = hi * "22 + ”23 * hi h5 * h 6 + ”27 + ”28 




(A 8 ) 
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If the equation for the downwash velocity at control point 1 is 
expanded as 6in exsurqjle, the result is 


^'^1 " ^11^1 ^12^2 ^15^5 ^15^5 ^16^6 




(A9) 


Alternatively, hy use of the relations in equation (A 5 ) , equation (A9) 
can be written as 


'^i " ^ 11^1 ^ 12^2 


^ 15^3 + ^4^4 ^84^5 ’’^ 85^6 


^ 82^7 ■'■ ^ 81^8 


(AlO) 


In case a symmetrical airload condition is being investigated, sub- 
stituting equations (a 6 ) into equation (A9) results in 


Wig = + ^2(^12 %?) ^5(^13 % 6 ) 


^ 4(^14 + 



(All) 


For an antisymmetrical airload distribution, substituting equations (AY) 
into equation (A9) results in 


^lA = ^1(^11 ■ ^18) + ^2(^2 “ % t ) ^5(^15 ■ ^16) 


F4(Ki4 - K 15 ) 


(A12) 
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^ (AI3) 


For the typical horseshoe system assumed, the complete series of 
equations relating downwash velocities to the circulations are: 

For a symmetrical airload distribution, 

''^Ig = ^1(^11 %s) ^2(^12 ^17) ^5(^15 %6) 

^2s = ^1(^1 + ^2S) + ^2(^2 + ^27) ^5(^23 ^26) + 

^4(^24 ^25) 

^3s " ^1(^51 So) + ^2(^32 + S?) ^3(^35 ^36) + 

^(^34 K55) 

''^ 4 g = ^i(^ 4 l ^ 48 ) ^2(^42 ^47) ^5(^43 ■*■ ^45) ■*■ 

F4(Kld4. + 

and for an antisymmetrical airload distribution, 

= h 4 ll - ''is) + ^s{\2 - ''17) + ‘'5(''l5 - >'16) 1 

r^CKii, - K15) 

"2;^ “ '‘ifel - ■'ss) + ^2(^22 ■ ''27) + - *'25) + 

Tljfei, - K25) 

"3a = '’l(''3i - >'38) + '’ 2(^52 - K,^) + r,(K 35 - Kjg) + 

- K55) 

"•lA ' ''l(''lH - ■'m) + ^ 2(^2 - ''47) + r,(K4j - K^g) + 

'’ 4 (K 44 - Kitj) 


>- (ai 4 ) 
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From equation (A1) , which expresses the relationship that must exist 
"between the downwash angle w/v at each control point, the wing angle 
of attack cxp, and the section lift-ciurve slope m^ for the wing station 

at the control point, the following series of equations result, where. 


If 


represents either 

Wt or w-i ; 

^A 


!5 = 

“5°3 

II 

l> 

V 

2 « 

^2 

W 4 _ 

“ 40-4 

II 

ro 

V 

2 jt 

^11 ~ 

[Kil t K^q) 

1 


^Ki2 t K^j) 

1 


(AI 5 ) 


or 


Hk = ^^(^44 i K 45 ) 


k = 4 «(Kl ± %) 


(AI 6 ) 


where the upper sign is used for symmetrical airloads, the lower sign 
is used for antisymmetrical airloads, and the subscripts L and R in 
equation (A16) mean left and right wing, respectively, then substitution 
of equation (A16) into equations (A14) gives 

— 

4 rtWj^ = + ^12^2 ^ 13^3 ^14^4 


4rtW2 

4rtWj 

4rtW2j^ 


^ 2 ± ^1 ^2 ^ 2 ^23 ^ 3 ^24 ^4 

^31^1 + ^32^2 + ^33^3 

Hl^l ^2^2 ^3^3 ^44^4 




-/ 


(AI 7 ) 
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The relation between the running load I and the circulation is 


I, = pvr. 


lo = 


It. = 


lu = 


PVF2 

pvr^ 

pvri^ 


(A18) 


and the equation relating dynEimic press-ure to mass density and true 
airspeed is 


q = pV^/2 


(AI9) 


Equations (AI5 ) , (AIT) , (AI8) , and (AI9) may be combined to give the 
following final system of linear equations: 


kiili + = l+qm^^o^ 


^21^1 ^22^2 ^25^5 ^24^4 ~ 4qm2Ct2 


+ ^j2^2 ^53^3 ^ 54^4 ~ 4qmj(Xj 


(A 20 ) 


hi’'! + ^2^2 H 5 S ^44^ " 
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Equations (A20) can now be written in matrix form. The equation 
for each station can be written as 


0=1 


kijlj = 4qm^o^ 


0=1 


^ ^5J^J ~ 

0=1 


0=1 


(A21) 


Equations (A21) and therefore equations (A20) can be expressed in the 
general form 



k. .Z . = 4qm.a. 
ij J 11 


(i = 1^ 2, 3, . . . n) 

(A22) 


In matrix notation equation (A22) becomes 

[s^ 0} = 

where 


ai 

OJ 

1 — 1 

^13 

ki4 

■21 

^22 

OJ 

OJ 

31 

^32 

^33 

^34 

■41 

ro 

H3 




(A23) 


(A2l) 
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and 



(A25) 


(A26) 


The development of the steps necessary to compute the elements 

in the [Sl] or "downwash" matrix follow. As a first step the relation 

which exists between the strength of a segment of a straight-line vortex 
and the induced velocity at a nearby point should be found. If the 
strength of the vortex (whose sense is given by the right-hand rule for 
moments) is V , the velocity induced at the point P can be written as 

r(cos a - cos b) 

wp = (A27) 

where a and p are the angles between the direction of the vortex seg- 
ment and lines joining the ends of the segment to the point as shown in 
the sketch 




Equation (A27) is used in the subsequent derivation of the 


matrix. 
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A plan view of the geometry of a typical horseshoe vortex on the 
left-hand wing is given in the following sketch, in which distances and 
angles are considered, positive as indicated and the sense of the circu- 
lation of the elements of the horseshoe vortex is given by the right- 
hand moment rule: 





The points V and D are a typical horseshoe reference point and a 
typical control point, respectively. 

The incremental downwash velocities induced by a single horseshoe 
vortex, if downwash velocities are considered as positive, are: 
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(1) For the left-hand trailing vortex the relations 


R = _(sy - h) 


a = 0^ cos a = 1 

P = 270^ - 0 cos p = -sin 0 

are substituted into the general relation 


_ r(cos g - cos p) 

to obtain the incremental downwash velocity 


r(l + sin 9) 

l.*(Sy - h) 


(2) From the right-hand trailing vortex where 


(A28) 


(A29) 


R = Sy + h 

a = 90° - 0 cos a = sin 0 

P = 180° cos 3 = -1 

the incremental downwash velocity is 

r(sin 0+1) 

^ 4«(Sy + h) 


(A3O) 
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(3) From the bound vortex where 


R = s. 


a = 180° - 0 cos a = -cos 9 


P = 180° - 0 cos 3 = -cos 0 


the incremental downwash velocity is 


Wb = 


r(cos 0 - cos 9) 
4ns„ 


The total downwash velocity at a typical control point due to 
complete single horseshoe is then 


W = Wj^ + Wj^ + Wr 


p / 1 + sin 0 1 + sin 0 cos 0 - cos 0 


J 4 -K \ Sy - h 


Sy + h 


Substitution of the identities 


Sy - h = 


Sy + h e 


Sx cos 0 
sin 0 

Sx cos 0 
sin 0 




into equation (A52) yields 


w = 


fl + sin 0 1 + sin 0' 


il-rts \ cos 0 cos 0 


(A51) 


(A32) 


(A 55 ) 


(A 54 ) 
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Equations (A 3 ^) and (A 2 ) are Identical; therefore. 


K 


ij 


1 /l + sin 
4 nsy. A cos . 


l_+_si^\ 

cos e^j ) 


If, as indicated in the sketch 


(A55) 



the control points are assumed to be located on the left semispan of 
the wing so that 6 l and ^ represent the pertinent angles for a 

horseshoe located on the left semlspan and and represent the 


pertinent angles for the corresponding horseshoe on the right semispan, 
then for a typical control point 




^ = 


1 /l + sin ^ 1 + sin 0j 


l4-rts^y cos ^ cos 0, 


Kr = 


1 /l + sin ^ 1 + sin 0j 


Ij-jts \ cos cos 0j 


(A56) 
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From equation (.AI 6 ) , 


k = t %) 


so that in the PO or downwash matrix each element k is coniputed 
from the equation 


1 fl + sin 0L ^ ®L ^ ^ 

k = — “ 


+ sin 6 


R' 


Sx\ cos 0L 


cos 0T 


cos 




cos 0 


U5T) 


R 


As in equation (A16) , the upper sign is to be used for symmetrical span- 
wise airload distributions, and the lower sign is to be used for anti- 
symmetrical spanwise airload distributions. Note that in equation (A34) 


in 0 = 


xp - Xy 


sin 


||(=^ - * (s'd - 5'v + *')“ 


COS 0 = 


yp - yy + h 


^(xp - x^)^ + (yp " yy + 


sin 0 = 


^ 




cos 0 = 


yp - yy - h 


\(‘=D - 
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The 


matrix therefore is computed from the matrix equation 


fSn 


1 p + Sin 0L 

1 + sin ©L 1 + sin ^ 1 + sin 0jA 

Li 


Sx\ cos ^ 

cos “ cos ^ cos J 


(A38) 


where the upper sign is used for symmetric flight conditions and the 
lower sign is used for antisymmetrical conditions. 

Since the r^ii matrix is used in equation (A 25 ) , which is 




the elements of the N matrix are seen to he influence coefficients 

relating the incremental downwash angle at each control point to the 
intensity of the running lift over each increment of the semispan of the 
wing. In general, all the elements in the principal diagonal of the 

matrix will always be positive and those elements not in the 

principal diagonal will always be negative because the velocities were 
considered as positive downward and wash velocities from a horseshoe 
vortex are downward only in the region behind the bovind vortex and 
between the trailing vortices of that horseshoe. 


Compressibility Corrections 

The method by which con5)ressibility effects are handled in this 
report is based on simple sweep theory. This theory is presented in 
references 6 to 8 and substantiated in references 9 to 12. Summarized 
briefly , infinite-aspect-ratio sweep theory indicates that compressi- 
bility effects are functions of the effective Mach mnnber Mg, which 
In this case is the Mach ninnber of the stream velocity component that 
is normal to the leading edge, so that the lift-curve slope in com- 
pressible flow is given by 


mo = 


m 

/l - M^cos^Am 


(A59) 
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where 


m 


section lift-c\irve slope at M = 0 


M 


stream Mach niimber 


Am 


effective sweep angle for compressibility effects or yaw 
angle of infinite wing, radians 


The same relationship exists for the effects of Mach number on 
section pitching -moment coefficient 


Compressibility effects on a finite-aspect-ratio swept wing can be 
handled in a similar fashion. For the finite-aspect-ratio swept wing, 
however, recognition should be given to the fact that compressibility 
effects will vary across the span of the wing. In general, smaller 
values of the effective sweep angle are indicated for the wing root and 
tip sections than for the mldsemlspan region. Even in the midsemlspan 
region the effective sweep angles for compressibility effects are not 
functions of the wing plan-form geometry alone; the spanwise variations 
of camber, thickness ratio, chordwise thickness distribution, and angle 
of attack are likewise involved in the determination of the value for 
the local effective sweep angles. 

Consider equation (7) for the downwash angle induced at a distance 
rearward of the lifting line which was derived from two-dimensional 
considerations 


This equation expresses the relation between the downwash angle at the 
three -quarter-chord point (i.e., at r = c/ 2) and the geometric angle 
of attack af when the section lift-curve slope is equal to its theo- 
retical value 2k . 

Several approaches can now be made to the compressibility problem 
regarding the values of lift-curve slope, the angles of attack, and the 
location of control points to be utilized for equating downwash angle 
to geometric angle of attack. In the method used in this paper the 
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distance r is always equal to one-half of the local streamwlse chord, 
so that in matrix notation equation ( 7 ) reduces to (see eq. (9)) 


r-0-1 



The applicable values of mo to be used in this equation are those 
at the effective Jfe.ch number Mg = M cos A^, and thus each wing station 
is permitted to have its own compressibility corrections . In the absence 
of test data reduced according to the method of appendix G, a value of 
mo can be calculated from equation (A59) for a value of Am equal to 
the geometric sweep angle of the plan-form strearawise quarter-chord line. 
Although, in general, this procedure will yield only an approximation 
to the correct theoretical value, this value will usually be sufficiently 
accurate for preliminary design purposes when A/cos A is large (where 
A is the aspect ratio) or M cos A is small, or both. 

This "modified angle of attack" method is characteristically dif- 
ferent from the "plan form distortion" methods of handling compressi- 
bility effects (refs. 13 and 1<4); the former distorts the angle at the 
original three-quarter-chord point, whereas the latter stretches the 
half-chord distsnice. 

An important advantage in the con^iressibillty-correction method 
presented in this paper is that only one [Sf] matrix need be computed, 
whereas a different (Sj matrix for each Mach number is required in 
the plan- form-distortion method. 

The following developments show the equivalence of these two methods 
of accounting for compressibility effects. The first comparison is for 
the case of an infinite wing in yaw. 

If Am Is the sweep or yaw angle for this infinite wing, then 
according to the plan- form-distortion procedure (indicated by the sub- 
script pd) the equivalent wing in incompressible flow (M = O) is to 
be at a yaw angle such that 


tan = 


tan Am 
/ l - # 


Further, the lift-curve slope of this equivalent wing in Incomipressible 

flow, multiplied by the factor 1//1 - M^, is equal to the lift-curve 
slope of the original wing in conipressible flow. 
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The lift equation for the equivalent Mach number method of this 
paper (indicated by the subscript is 


^Me " 


2n cos Am 
^ - W^cos^M 


and that according to the plan-form-distortion theory is 


Lpd = qca 


2jt cos Am 
/l - 


In order for the two lifts to be equal, the following equality must be 
true for any value of M: 


cos Am cos Am 

l/l - M^ ^1 - M^cos^Am 

It can be shown that such is the case by the use of elementary trigo- 
nometric identities. 

The next coniparison is for the case with zero sweepback. The plan- 
form-distortion method stretches the half-chord distance as a function 
of Mach number to 

I -1 

£— = 2. 1 

which when substituted for r in equation (3) gives 


ry/i - 

JtC 




W 
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Equating the downwash angle at the stretched control point to the 
section angle of attack gives 


w 

V 


ri/l - 
itcV 


= Of 


The compressible lift-curve slope (ref. l4) is 


(a4i) 


mo “ 


2rt 

]/l - 


(a42) 


Equations (a4i) and (a 42) may be substituted into equation (2) to give 

2n r/l - m 2 pv^c 
|/l - 2 

or 


I = pvr 

For the compressibility-effects method of this paper the downwash 
angle at the three-quarter-chord point equated to the angle of attack 
corrected for effects of section lift-curve slope is (see eq. (8)) 


[v\ _ r _ ^ 

VV/3C/4 «cV 2it 


from which 


Of = 


2 F 

moCV 


(a43) 


When equation (a 43) is substituted into equation (2) the same final 
answer is obtained; 


1 = pvr 

The identities of compressibility corrections in the case of either 
the unyawed (Aj^j = O) infinite wing or the yawed (Aj»j O) infinite wing 
substantiates, at least qualitatively, the use of the compressibility 
correction methods as presented in this paper. 
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APPENDIX B 


THE ELASTICITY MATRICES QSgJ AND [p2'l 


The purpose of this appendix is to explain, in brief form, the 
structural theory involved in and the steps necessary for the con^juta- 

tion of the elasticity matrices and 


Development of the 



Matrix 


In the development of the downwash matrix in appendix A, a 

continuously varying spanwise airload distribution was replaced by a 
series of constant-intensity rimning loads, each of which covers an 
increment of the wing semispan. This concept of an equivalent system 
of loads is also used in this appendix. 

Consider the geometry of the structural skeleton of the wing to be 
as represented in figure 5, in which double arrows indicate that the 
right-hand rule of moments applies. The notation in figure 5 is defined 
as follows: 





d^, dg, . • • 


total lift of increment of wing having span of 2 h^^, 
2h2, . . • 2 hj^, respectively, numbered inboard 
from the left wing tip, lb (see eq. (Bl)) 

strearawise distance from horseshoe reference point 
at a wing station to elastic axis at the same 
station, positive when elastic-axis point is to 
rear of horseshoe reference point, in. 

strearawise distance from elastic-axis reference 
point at a station to elastic-axis reference 
point at the next inboard station, positive when 
inboard elastic-axis point is forward of outboard 
elastic-axis point 

lateral distance between elastic-axis point at a 
station and elastic-axis point at the next station 
inboard 

local sweepback angle of elastic axis 
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%2’ 

%2^ 



rolling moment at elastic-axis point around longi- 
tudinal axis through local elastic-axis reference 
point due to total lift of all the vortices out- 
board of this point, positive when it raises left 
wing tip, in- lb 


. My pitching moment at elastic-axis point around lateral 
axis through local elastic-axis reference point 
due to total lift of all the vortices outboard of 
this point, positive when nose up, in- lb 


M, , Mp, . . . beam bending moment at elastic-axis point about an 

axis perpendicular to local elastic axis, positive 
when it puts compression in wing upper surface, 
in- lb (see eq. (B2)) 


Ti, Tp, . . . Tjj torsional moment axound elastic axis at local 

elastic-axis point, positive when it is in direc- 
tion of leading edge up, in-lb (see eq. (B3)) 


The general form for the wing lift is 


K (Bl) 

Where is the intensity of the running lift at station measia'ed 

in po\mds per inch and 2h^ is the span of the horseshoe vortex at 
station the general form for the bending moment is 


cos An - sin An 


(B2) 


and the general form for the torsional moment is 




n 


cos An + M^n 




(B3) 
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At station 1, on the center line of the horseshoe vortex nearest the 
left wing tip (see figs. 3 and 4) , the following equations apply: 




Li hj^ Lj^hj^ 


(Bl^) 


Mil =Tfl 


L]_ / hj^ tan A]_ 


) 


(B5) 


which, when substituted into equations (B2) and (B3) , yield 


Ml = L^l 


hj^ ej^ sin 


(4 cos A^^ 


(b6) 


Ti . Lj, 


cos 


(B7) 


At station 2, on the center line of the second horseshoe vortex 
inboard of the left wing tip. 


%2 ~ 


+ 


^2^ 

T” 


(B8) 


%2 "" ^l(®l " 



"tcin A^ 


(B9) 


H 
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and equations (B2) and (B 3 ) become 


^2 = cos A 2 - - fi)sin Ag + Lg^- 


hg 62 sin Ag^ 


^ 4 - cos Ag 2 


(BIO) 


Tg = Lj^dj^ sin 


Ag + 


^ cos Ag "i" Bg 


eg cos Ag' 


(Bll) 


At station 5^ 


% =■ h{h + '^ 2 ) + V2 


Ljhj 


(B 12 ) 


%3 = hh - ‘■1 - ^2) * ^2(^2 - ^2) * - 


tan A^' 


and, similarly. 


”5 - + vj ] ^ - fi(®i - fi - ^ 2 ) 


(BI 3 ) 


Lg ( e 2 - fg ^ 


6 o - fo) sin Ai + L 


sin 


3 ^31 


cos 




(Bl4) 


Tj — t Lgd^sin ~ fq ~ ^"2) 


(®2 " 


cos Aj + Lj 


cos A^ 


( BI5 ) 
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Equations similar to these can be developed for M and T for 
the remaining stations on the left semispan of the wing. 

In algebraic form, the equations for and Tj^ are 



L 


n 


/ 

W cos 


e^ sin An' 


(BI6) 



Ln 


^n ^ 
2 


(BIT) 


Notice that the first tvo terms in each of these eqiiations are equal 
to zero vhen n = that is^ for the tip station. 

In matrix notation, equations (BI6) and (BI7) "become 


{«}- 

{t} = 


E°® ElM - ^][“] 




in JrgJ + [cos [V] -[hj- 


(BI8) 


(BI9) 
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where 


R = 


Ml 

M, 


Mn 


{T} = 


w- 


T, 


■-n 


Ln 


(B20) 


(B21) 


(B22) 


[cos = 


cos 

0 

0 

0 


[sin ^ = 


sin Ai 
0 
0 
0 


0 

cos A2 
0 
0 


0 

sin Ag 
0 
0 


0 

0 

cos A5 
0 


0 

0 

sin A^ 


0 

0 

0 

cos Aj|^ 


0 

0 

0 

sin Aj^. 


(B25) 


(B24) 


0 


0 


0 


0 0 



*01 

o 


(B25) 


it- cos^A^ 


[fi] = 


d-i + 


+ d^ 


di + dg + d^ + d|^ 


4 COS^Ag ^ 

dg + d^ 

d 2 + d^ + d|^ 


h^ e^ tan A^ 


4 COS^Att ^ 


d^ + di^^ 


h4 tan A^ 


4 cos^Aj^ 2 


h^ e^ tan A^ 

4 cos^A^ 2 


(B26) 
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P2] - 


2 tan 
<^1 

^2 ■^" + ^Ij. 


2 tan ^2 
d2 

d-2 ^3 

d2 + dj + d]^ 


0 


2 tan A^ 
^3 + <3-i4. 


0 


0 


e^ 


2 tan Aj^ 
d4 


0 


2 tan 




(B2T) 


Equations (B18) and (B19) provide the means for confuting the Bending and torsional moments 
along the span of the wing due to the loads L as given by equation (B1) . 

Each of the loads I or L, however, is affected by the variation in angle of attack from 
root to tip. This problem is handled by computing the streaniwise angle -of -attack change ag 

due to the structural deformations caused by the series of loads {l}. Since a strearawise angle- 
of -attack change is required, assume that a unit positive pitching moment m^ (nose-up moment) 

is applied in turn at each of the points 1, 2, ... n on the elastic axis. The unit pitching 
moment is in the plane in which Og is to be measured and its direction agrees with the posi- 
tive sense of ag. 
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The streamwise angle -of -at tack change can then be obtained by- 
applying the general relation (see ref. 15 ) 


aso = 


0 


mM ds 


tT ds 
GJ 


(B28) 


where 

ag desired streamwise angle -of -attack change due to values of 

^ bending moments M and torsional moments T along elastic 

axis of wing resulting from the series of loads -{^L} 

m beam bending moment per unit pitching moment applied at sta- 

tion at which as^j is to be determined 

t torsional moment around elastic axis per unit pitching 

moment 

ds increment along elastic axis 

El effective beam bending stiffness around the axis of the bending 

moments M and m 

GJ effective torsional stiffness around the axis of the torsional 

moments T and t 

The stations on the wing for which the angle-of -attack changes ag 
are to be computed are those on the center line of each horseshoe vortex. 

It is assumed that sufficient accuracy in the results will be 
obtained by using the values of bending moment M and torsional 
moment T obtained by means of equations (Bl8) and (B19) and that these 
values of M and T may be considered to be constant over the incre- 
ments of the wing span corresponding to the span (2h) of each horseshoe 
vortex. 

Values of El and GJ are also assumed to be constant over each 
svich increment in wing span, and the values to be used are the effective 
values which correspond to the wing station t) at the midpoint of each 
horseshoe vortex. 

The general relation (eq. (B28)) can be used to express ag in 
matrix form by letting 
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Mn 


sweepback angle of elastic axis at station radians 

beam bending moment at station r|j^ due to loads on wing 


n 


torsional moment around elastic axis at station due to 

loads {^li} on wing 






beam bending moment at station i per unit pitching moment 
applied at station j 

torsional moment around elastic axis at station i per unit 
pitching moment applied at station j 


(ds)^ 

2hn 


2h 


h 


COS A, 




horseshoe span at station t]. 


'n 




'n 


angle -of -attack change at station due to all the 

loads on wing 


Then, for station 1, 


a„ = 

®1 




(EI)-j^ cos Aj_ (EI)2 cos A2 




(EI)^ cos Aj^ 


11 ^ 1^1 


+ + 


(GJ)]_ cos (GJ)2 cos A2 


^^nl'^n^ 


(Cy)j^ cos An 


(B29) 


for station 2, 


a, 


S2 


^2^l\ 


m^^Mghg 
+ + 


(EI)-]_ cos Aj^ (EI)2 cos A2 




(El)n cos An 


2ti2Tihi t22^2^ 

+ + 


cos Aq^ (01)2 cos A2 


^^n2^n^ 


(GJ)n cos An 


(B30) 


and so forth. 
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From the geometry of the problem, when point i is at or inboard 
of point j, the bending moment due to a unit pitching moment applied at 
point j is 


mij = -sin Ai 


(B 31 a) 


when point i is outboard of point j, the bending moment is 


m. 




= 0 


(B31b) 


when point i is at or inboard of point the torsional moment due to 
a unit pitching moment applied at point j is 



cos 


(B32a) 


and when point i is outboard of point j, the torsional moment is 



(B32b) 


and the complete series of equations for ag, of which equations (B29) 
and (B 30 ) are representative examples, can be written in matrix form as 



o 

~2h 
cos A 



o 



2h 

cos A 



(B33) 
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Typical elements in [mj and Qt^ are 



-sin Ap/2 

-sin A2 

-sin A, 
3 

-sin Aj^ 


0 

-sin A2/2 

-sin A^ 

-sin Ai^. 

II 

0 

0 

-sin A^/2 

-sin A|^ 

0 

0 

0 

-sin Aj^y'2 


(B5U) 




cos Aj^/2 
0 
0 
0 


cos A2 
cos A2/2 
0 
0 


cos A^ 

cos A^ 

cos 

cos Ai^. 

cos Aj/2 

cos Aj^ 

0 

cos Aj^/2 


(B35) 


Substituting equations (B1) , (BI8 ) , and (BI9) into equation (B33) 

gives 




H 


o 

2h 

cos. A 


0_- r- 
1 

El 


[cos ^Jfl] - fin^[u] 


o o ^ 




2h 

cos A 


1 

GJ 


V-» KJ 

l^in + [cos 


|2°h]{z} (B36) 
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If 


II 

0 

r-if 

[mj 

-^[_cos A 

0 

r 1 

[ ei _ 

0 0 

[cos A^[rJ - Qsin ^[1^ 





0 

"“ll 

GjJ 

~ 0 0 ^ 

0 in ^ -^[uj 


(B37) 


then 


h) - MO} 


(B38) 


Eq.uation (B5T) represents the most general form for the matrix, 

and each element a . . of this matrix represents the angle-of -attack change 

in radians at station i due to the structxrral deflection of the wing 

caused by a unit loading at station j. In effect, the matrix is 

an array of influence coefficients, and the elements of this matrix may 
be con^juted according to equation (BJT), or, when an actual wing is avail- 
able, they may be obtained by load-deflection tests of that wing. 


Development of the Auxiliary Elasticity Matrix 

The matrix as defined in the previous section considers air- 

loads outboard of each specified station, as defined by a strearawise cut 
through the wing box. Inasmuch as ribs are conventionally installed in 
a plane essentially normal to the elastic axis, the spar -box loads are 
more correctly determined by considering the airloads outboeird of a plane 
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nonaal to the elastic axis. These loadings are represented in the fol- 
lowing sketch: 



Elastic 

axis 




The corrective loading is then the loading obtained by subtracting 
the loading of @ from that of (S) . The trianguleir areas (c) aind (5) 
represent the corrective loading. It will be noticed that the loading 
over the area (c) is the negative of the loading of the corresponding 
area in @ and the loading over is the positive of the corresponding 
area in @ . In order to obtain these corrective loadings^ it will he 
necessary to assume a distribution of pressure over the corrective areas. 


The equations are now developed that represent the shear, wing 
moment, and wing torsion due to the corrective loadings, and the 

matrix is modified to include these effects. A plan view of the 
left wing is shown in figure 5 ^ where 


E local angle between elastic axis and lateral axis equal 

to A, the sweep angle of the elastic axis 

L local angle between leading edge and lateral axis 

T local angle between trailing edge and lateral axis 


The correction terms for moment, shear, and 


airload in the 


matrix are 


torsion due to wing 
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• • • 

• • • ^n 

m^, ^2, 

• • • ^ 

ATi, 


. /^. 


n 




correction to rolling moment at elastic-axis point 
around longitudinal axis throiagh local elastic -axis 
reference pointy positive when it raises left wing 
tip^ in- lb 

correction to pitching moment at elastic-axis point 
around lateral axis through local elastic-axis 
reference point, positive when nose up, in-lb 

correction to wing bending moment at elastic-axis point 
about an axis perpendicular to local elastic axis, 
positive when it puts compression in wing upper 
surface, in-lb (see eq. (B59)) 

correction to torsional moment around elastic axis at 
local elastic-axis point, positive when it is in direc- 
tion of leading edge up, in-lb (see eq, (b4-0)) 

correction to local shear, positive up, lb 


. ^ 


n 


The wing bending moment is 

Mn = cos En - sin 


and the wing torsional moment is 


^n = Ej„ + sin E^ 


(B39) 


(b40) 


Assmed pressure distribution .- The distribution of pressure p 

over the corrective areas is assumed to be given by the two-dimensional 
equation for a thin flat plate 


p = 4qa' 


1 - 

c 


x/c 


(b4i) 


and modified to account for the actual value of local running lift. In 
the modified form^ the distribution of pressure (measured in poimds per 
square inch) along a streamwise chord becomes 


2 I 

p 


n c ^ 


1 - 


x/c 


(Bl)-2) 
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where x/c is the fraction of the strearawise chord to the rear of the 
leading edge. The ratio z/c is assumed to be constant so that the lines 
of constant pressure radiate from the origin 0 in figure 5. In accord- 
ance with the assvmrptions previovisly given concerning the corrective 
loadings, the pressures are negative (^downward) over area FDG and posi- 
tive (upward) over area ABG. 

A polar coordinate system (p,0) is selected, and the shear, moment, 
and torsion due to the pressures on area OFB will be subtracted from 
the shear, moment, and torsion due to the pressures on area ODA. Lines 
of 0 = Constant are lines of constant pressiire and p will be a func- 
tion of 0 only. 

From the geometry of figure 5, p may be expressed as a function 
of 0 by writing 


Hence 


and 


or 


c + OM tan T - x = OM tan(L - 0) 


X = c + OMjtan T - tan(L - 0)^ 


c + OM tan T = OM tan L 


(b43) 


OM = 


tan L - tan T 


(m) 


Combining egimtions (b 43) and yields 


X = c 


1 + 


tan T - tan(L - 0) 
tan L - tan T 


(b45) 


If equation (b 45) is substituted into eqimtion (b42), 


p = K 


tan(L - 0) - tan T 


tan L - tan(L - 0) 


(b1^6) 


K = 2I 
nc 


(bi+t) 


where 
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Rolling-moment correction at station n, The rolling-moment 

correction at station n, is obtained by adding the rolling 

moment due to area OFB to the rolling moment due to 

area ODA (fig- 5)* 

The rolling moment due to area OFB may be found by writing 

the equation for the elementary moment about FM for any value of p 
between 0 and that at the line FM as : 




dAMjf = pj^M - p cos(L - 9^ p dp d0 


The eqiiation of the line FM is 

p cos(L - 9) = OM 


or the value of p at the line FM is 


OM 

cos(L - 9) 


Over area OFB the moment is negative; hence, 


— -KOM 


*L-T 


0 


tan(L -9) - tan T 


tan L - tan(L - 9) 


d9 


OM 


cos 


(L-9) 


P dp + 


0 


OM 


K 


PL-T 


0 


tan(L - 9) - tan_T _ e)d9 

tan L - tan(L - 9) 


cos 


(L-e) 


p^dp 


0 


Since 9 varies between 0 and L - T, the integrals can be evaluated 
by letting 


tan L - tan(L - 9) _ 
tan L - tan T 

t an(L - 9) - tan T _ po-t-^n 
tan L - tanCL - 9) “ 




(b48) 
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It follows that 


d(L - 0) = -2(tan L - tan T)cos^(L - 0)sin U cos U dU 
Consequently^ 


= -KOM^(tan L - tan T) cos^U dU + 


^ OM^(tan L - tan T) T ^ 

5 Jo 


cos% dU 


= -IlLKOM^(tan L - tan T) 
12 


lc‘ 


6(tan L - tan T)^ 

The rolling moment area ODA may be derived in 




similar manner to get 
PL-T 

= KOM 




tan(L - 0) - tan T 


tan L - tan(L - 0) 


d0 


cos(L-0-E) 


p dp - 


0 


PL-T 


K 


UO 


where 


cos(L - e)ae 

y tan L - tan(L - 0) 


cos(L-0-E) o 

P^dp 


(B50) 


_ _ OM cos(E - T) , , 

R - + c(l - e)sin E 


(B51) 


cos T 
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After integration this equation simplified to 


rtK0MR2(tan L - tan T) 

" 

4(tan L + cot E)^sin^Eljl - 


nKE5(tan L - tan T)(4 - 3Q^) 

■ (B32) 

24(tan L + cot E)^sin5Eu(l - q2)3 


where 


2 tan L - tan T / ^ P \ 

= (0 < < 1 ) 

tan L + cot E 


Equations (BJf9) and (B52) may now be combined to give the rolling- 
moment correction at station n: 


^ A 


•B 


jtK 

6cR^(l - Q^)(tan L + cot E)sin E - R?(l- - 3Q^)(tan L - tan T) 

2k 


j(l - Q^)?(tan L + cot E)?sin?E 


2c? 

(tan L - tan T)^ 


If U defined in equations (bJ+ 8) is written as 


cos(e - T)cos L 
sin E sin(L - T) 


(B55) 


(B5U) 
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and if 


V = tan L + cot E 
then R (see eq. (B^l) ) becomes 

R = c(U - e)sin E 


(B55) 


(B56) 


and 


^Xn 


Ic^ 

6(U - e)^(l - q 2 ) . (u - e) 5(4 - 5Q^)q2 

12 

V^(l - 



(B5T) 


For the special case of the untapered wing where the leading edge is 
parallel to the trailing edge or T = L, 


l£ sin^E c o s ^ + 1) 

16 cos^(E - L) 


(B58) 


Pitching-moment correction at station n, AMy^.- The incremental 

pitching-moment corrections due to the areas OFB and ODA may be 
derived and combined in a manner similar to that for the rolling-moment 
corrections and the equation for becomes 

q 2 (u _ e )5 1^1,. L - + 6 ( 1 - Q^)(U- e)^(€Q 2 v - tan L) 

V2(i _ q2)5/2 


ic^ 

= -< 

^ 12 


(6e -l)VQ^- 2 tan L 
v2q4 


(B59) 


For the special case of the untapered wing where the leading edge is 
parallel to the trailing edge or T = L, 



_ Iff (8e2 - lv€ + 1) 

16 


sin E cos L 
cos(E - L) 


2 


sin E cos L tan L 


(b6o) 


cos(E - L) 
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Shear correction at station n, ASn*- The elementary shear is 


dASn = PP dp d 9 


(B61) 


From this equation the incremental shears over the areas OFB and 
ODA may also be derived and combined in a manner similar to that for 
the rolling moment to give 


■"n 


Is. 

q2(u - e)2 

1 

2 

V^l - 



(b62) 


For the special case of the untapered wing where the leading edge is 
parallel to the trailing edge or T = 


.c. sin E cos L /, ^ x 

ASn = - ^ : r- - l) 


4 cos(E - L) 


(b 63 ) 


Modification of [33] matrix.- The matrix is modified by 

using the corrective rolling and pitching equations. Using the equations 


^ cos En - sin E„ 


(b 64 ) 




and letting (from eqs. (B57) and (B 59 )) 

P = -^n/ 
J = 


I 


(B65) 


(b66 ) 

(B67) 


1 


NACA TN 5050 


65 


results in 


= ("M- cos + J sin Ej^^) I 


ATn = -(^ sinEj^ + J cos Ejy)l 


or, in matrix form (with = Ej ^) , 


= 


00 00 

jcos + [sin 


0 


o o r* 1 

AT|- = - |sin a]H . [cos aJM W 


o 


Eqiiation (B18) thus becomes 


{“]- 


[cos a] [r[] [2h] - [w] 


- f a] H&'O - P] 




and equation (BI 9 ) becomes 


(t} = 


[s in ^ |t2^ ” M ^ ^ 


[u][2^ - [j] 




The corrected form of the 




matrix thus becomes 




M 


o 

2h 


cos A 


r_o 
1 


El 


o 

[^os 


[r[]ph] - H 


[sin 


[uj[2hj 


' — ' c 


2h 


cos A 


GJ 


[sin A^ 


f 2] P”! - p] 


(cos 


[M][2hJ 


(b68) 

(b69) 

(BTO) 

(B 71 ) 

o”l 

- p]j 

o"l 

-hJJ 

(B 72 ) 
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where, as in equation (B58) , 




(BT 5 ) 


It should be noted that, because of the assumed pressure distribu- 
tion given by equation (b 4 i), this correction does not correctly handle 
conditions involving deflected flaps and roll-control devices. 
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APPENDIX C 

COMPUTATION OF MATRICES 

This appendix outlines the method by which the component parts of 
the |a^ matrices may be determined. The matrices comprise two 

classes of twist: (l) those which would be present even if the wing were 

rigid and (ll) those due to inertia effects, thrust or drag, and section 
pitching moments on the flexible wing. A given airplane design condition 
may obviously require any combination of the twists listed under classes I 
and II, that is, 

{ag} = [agj^ + |ag|^^ (Cl) 

For simplicity, the sources of summarized as follows: 

Class I - aerodynamic twists (l.e., zero-lift-line shifts due to 
effects other than wing structural deflections) 

(a) Built-in twist due to camber or construction, or both. 

(b) Interference twist due to fiiselage, external stores, and 
so forth. 

(c) Twist due to flap deflection. 

(d) Twist due to aileron deflection. 

(e) Twiat due to spoiler deflection. 

(f) Apparent twist due to airplane rolling velocity. 

(g) Apparent twist due to airplane pitching velocity. Angles of 
attack due to airplane pitching velocity should be measured at 3c A. 

Class II - structural twists due to wing deflections caused by the 
following inertia and aerodynamic loadings which are independent of the 
wing lift distribution; 

(a) Vertical acceleration acting upon dry- wing dead weight, 
wing internal- fuel dead weight, and external- stores dead weight. 

(b) Effect of airplane rolling acceleration upon dry-wing dead 
weight, wing internal-fuel dead weight, and external- stores dead 
weight . 
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(c) Effect of airplane pitching acceleration upon dry -wing 
dead weight, wing internal-fuel dead weight, and external-stores 
dead weight . 


(d) External- stores thrust or drag. 

(e) Section pitching-moment coefficient with flaps, ailerons, 
spoilers, and so forth in neutral position. 

(f) Incremental section pitching -moment coefficient due to flap 
deflection, aileron deflection, or spoiler deflection, or in various 
combinations . 


Of the class I twists, the aerodynamic built-in twist is known from 
the wing geometry and the characteristics of the profiles used, and the 
interference twist may be obtained from actual tests of the complete air- 
plane configuration, either model or full-scale, by using the procedure 
of appendix G or from other calculation methods . Twists due to control 
deflection are determinable from the same type of tests as were used to 
0Y3,3.uate the interference twist . Apparent twists due to airplane rolling 
and pitching velocities are completely determined when the airplane flight 
condition to be investigated is specified. 


The type of twists due to the effects of wing deflections arising 
from loads which are independent of wing angle of attack, such as those 
listed under class II, may be computed with the aid of eqixation (B35) 


since, in this case. 






(C2) 


where -fw and are the wing bending moments and the torsions along 
the wing elastic axis due to the loadings of class II or to any combina- 
tion of them. 


Although the {"a 1 twists do not all explicitly require the cal- 

L 

culation of loadings, they do influence the equilibrium airload and must 
be accounted for in determining the net wing loads from equations (21), 
(22) , and (25) • 
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APPENDIX D 

DERIVATION OF EXTERNAL-STORE MATRICES 


The term external stores as used in this report is intended to apply 
to such items as nacelles , external fuel tanks, bombs, rockets, and 
similar items commonly attached to the wings of airplanes . > 

In this appendix the lift coefficients and pitching-moment coeffi- 
cients for each external store in the presence of the rest of the air- 
plane configuration are assumed to be known and to vary linearly with 
angle of attack. On this basis a set of linear equations is developed 
in which the airloads on each external store are accounted for in com- 
puting the deflections of and airloads on the elastic wing or airplane . 
The coefficients used are based on the airplane wing area and the wing 
mean aerodynamic chord. There is a specific lateral reference axis for 
each external- store pitching -moment coefficient, and the angle of attack 
of the external store is taken to be the angle of attack of that wing 
station within the span of whose horseshoe vortex the station of the 
external store lies . Special symbols used in this appendix eire as 
follows : 


pitching-moment coefficients, respectively, of an 
external store measured at its aerodynamic center (see 
eqs. (Dl)) 




lift and pitching-moment coefficients, respectively, of an 
external store at zero angle of attack measured at its 
aerodynamic center 


^ 0 -’ ^0 






% 


% 


e+1' 

%n 

e+l' 

%n 


lift and pitching moment, respectively, of an external store 
meas^lred at its aerodynamic center (see eqs. (D2)) 

rolling moment at elastic-axis point aroxind longitudinal 
axis through local elastic-axis point due to loads on a 
particular store outboard of this point, positive when 
it tends to raise left wing tip, in-lb 

pitching moment at elastic-axis reference point around 
lateral axis through this reference point and due to 
loads on a particular store outboard of this point, 
positive when it tends to raise nose, in- lb 


Additional symbols required are shown in the following sketch; 


TO 
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Consider one external store whose characteristics are 


so that 






'0 


(Cn^)^afE 


%0 = <lS jcLg^ + f LE)a°^^ 
^%:o = <lSc[Cj^^ + 


(Dl) 


(D2) 
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The lift and moment on the external store L^q and M^q 

with respect to point e hy the dimensions y^ and Xg 
the sketch. 

At the first station inboard of the external store 


Mg = cos Ag - My^ sin Ag 

= Leo (ye sin Ag) - sin Ag 


le = %g + %g cos Ag 

= Leo ( ye sin Ag + Xg cos Ag) + N^o cos Ag 
At the next inboard station e + 1, 

= LEo(ye + dg) 

^e+1 ~ ^ 1^0 (^€ " ^e) 

Me+1 = Leo [^ye + dg)cos Ag+p - (xg - fg)sin Ag+i^ - Mgo 
^6+1 = %:0 + <ie)sln Ag+p + (x^ - fg)cos A^^.^ + MgQ 


are located 
as shown in 


(D3) 

(Div) 


sin Ag+p (D5) 
cos Ag^o^ (d6) 
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At station e + 2, 

%e+2 "" I^E0(ye + <ie + <ie+l) 

^e+2 " ^0 Leo(x£ - fe - fe+l) 

Mg+2 = Leo [[ye + tie + <ie+l)cos Ag^.2 - (x^ - fg - fe+i)sin - 

%0 sin Ae^_2 (DT) 

Te+2 " Leo l^e + tie + ti£+E)sin Ag+2 + (^e ~ i"e ~ fe+l)cos Ae+^ + 

MeO cos Ae+2 (d8) 


At station the most inboard station considered, 


Since 



Mn = cos An - 
Tn = %n + %n 


sin Aj^ 


^n 


cos 
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then, 


Mn = Leo 


and 


'■ ^EO 


cos An^e - r dm^ - sin An(xe - 


sin An(ye + 


n-1 

5 : 

m=€ 


dm) + cos An(Xg 


n-1 


m=e 


'•m 


- %o sin An 


(D9) 


n-1 

m=e 


‘-m 


+ Mqq cos An 


(DIO) 


In matrix notation, the following equations, typical for each 
external store on the semispan, are equivalent to the algebraic equa- 
tions (D9) and (DlO): 


{%} 

(te) 



r 

[sin 


r ^ 

[cos 


'{l - 

(u } I ^EO + 


5in aJ [ie} %o (Dll) 
[^os a] [ie] Meo (D12) 


where, if the value at the top of the column is for the tip and the 
other values are for the successive stations inboard to the root and 
e is the station within whose horseshoe span the external store is 
located. 


0 

0 

0 



< 


ye 


ye 

ye + dg 

y^ + dg + dg4-i 
+ de + dg+i + dg^2 




(DI 5 ) 
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T4 



0 

0 

0 



Xe - fe - fe+1 
Xe - fe - fe+1 ~ fe+2 





(DI5) 



0 

0 

0 


; Mg V 
Mg +1 
Mg +2 


J 


(DI6) 



r 

0 

0 

0 


Te+l 

Te+2 


(bit) 


From equation (B35) 




r o ■ 
2h 

cos A 



+ 



(di8) 
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where the external-store bending moment and torsion are given by equa- 
tions (Dll) and (D12), respectively. 


to; 


Substitution of equations (D2) into equations (Dll) and (D12) leads 



Eos A]{lE}qSc(C„5.)^afg 


For convenience, let 


{^2m} = 
{^2t} = 




then equations (DI 9 .) and (D20) become 


(D19) 


(D20) 

(D21) 

(D22) 



(D23) 
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{Te} - ^]{%}‘=Ci"Eo * 

[cos a] {lE}=(Cn5)^OfJ 

Now let 

{kim} = {x2m}cLeo ■ ^ {^Ej^CmEo 

[k2m} = [x2M}(CLE)a ' ^ {^e]c (Cmg)^ 


|Kit} - |x2t|^Leq + |f°® A]{^E}cCmE^ 


{K2t| = (^2T}(CLE)a + a] {^Ejc (Cn^)^ 


then 


{Me} = qs|[KiM} + [K2M]afEj 

+ {^arjo^E j' 

Equation (D18) relating with 


c 


written as 


{«s}e - 1s||ki| + |K2ja^ 


E 


(D24) 

(D25) 

(D26) 

(D27) 

(D28) 

(D29) 
(D30) 
now he 

(D31) 
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where 


{Kx] 



O O 


2h 

cos A 



o o 




(D32) 


The equations just derived are typical for a single external store 
on the semispan. Each additional external store on the semispan requires 
a similar set of equations. 

The total angle -of -attack change due to the wing airloads shear and 
the airload on the external store is 


M ■ (Mw * {“=]e 


(D55) 


where the subscript W denotes wing and E denotes external store . 

Here there is an 
nal store; that is, 


He 


term on the right-hand side for each exter- 


N ■ N„ ^ ^ {“=1e2 * • • • 


(D54) 


Equations (l2), (l5)> and (B38) for 
repeated for convenience are 





respectively. 


78 


MCA TN 5050 


Including the effects of external stores on according to 

equations (D5^) yields 


o 


1 

4qmo 


[^1 (^] - {^ r } + + 



(D55) 


From equation (D51) 


{Mb, = 

We 2 ' ■" {‘^Is'^^Eg 


(D56) 


However, cLf can be found from. 
E 


“^En ■ LM{M 


(D57) 


where a typical value for given hy 



0 1 0 0 0 0 



(D38) 


The row matrix j_EnJ many colianns as there are horseshoe 

vortices on the semispan and has a value of unity entered in the column 
which corresponds to that horseshoe within whose span the external store 
is located, (in the example implied in eq. (D38) there are ten stations 
and the external store n is within the span of the fourth horseshoe 
from the left wing tip.) 

Combining equations (l2) and (D37) results in 




o 




(D39) 


NACA TN 3050 


79 


Combining equations (D35), (D^ 6 ), (D39), and (B38) yields 


0 

1 


4qmo 


- (“r} + {o-g} + [S2]{4 + 


qS 


KIiN 

0 

~ 1 ~ 
4qmo 

[Sl] {*}} 


0 

r 1 ~ 

4qmo 

[si]W' 


(DljO) 


Let 


[An] = {k2],L\I {%} = K)„ 


(D4l) 


and 




(Dll- 2 ) 


If eqxiations (d 41) and (d 42) are substituted into equation (d 40), 
the following matrix equation is obtained for ein elastic wing having a 
number of external stores on the semispan: 




[l} = {oj.} + [ag} + qS{B] 


(D43) 


This equation is similar in importance and usefulness to equa- 
tion (16) except that the aeroelastic effects of the airloads upon the 
external stores are now included. It reduces to equation (16) if all 

o 

the I^A^ and matrices are set equal to zero, since Til = [j^'^ 
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APPENDIX E 

WING-FUSELAGE INTERFERENCE 

In this appendix an approximate method of calculating the effect 
of a fuselage on the spanwise airload distribution on the wing is devel- 
oped. The primary so\arces of this interference are considered (ref. 16) 
to be as follows; 

(1) The effects of the image vortex system inside the fuselage on 
the downwash angles at the various control points over the wing semispan. 
This image vortex system induces a flow which is a first approximation 
to that necessary to satisfy the condition that there be zero velocity 
normal to the fuselage . 

(2) The effects of the vertical "overvelocities" resulting from the 
fuselage angle of attack. These velocity increments affect the local 
angle of attack at the various control points on the wing. 

This method of analysis can be used when applicable data of the type 
described in appendix G are not available . 

Item (1) may be considered as the effects of wing airloads upon 
themselves due to the presence of the fuselage; whereas item (2) accounts 
for the effects resulting from the fuselage having an angle’ of attack. 

In the development that follows the fuselage is ass'umed to be of 
circular cross section, of constant diameter, and infinitely long. 

For the image vortex system mentioned in item (l) , the individual 
images of the wing trailing vortices can be sho'wn to be located on a 
straight line joining the axis of the fuselage with the axis of the par- 
ticular wing trailing vortex at a distance from the fuselage center line 
such that 


R 


2 - 



(El) 


where a is the fuselage radius and R^^ is the distance from the fuse- 
lage axis to the trailing vortex (see fig. 6) . 

In similar fashion the bound vortex is assumed to have an image 
within the fuselage cross section. The image of the bound vortex is 
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assumed to lie on a straight line joining the forward ends of the image 
trailing vortices . The forward ends of these image trailing vortices 
are assumed to lie in the same transverse plane (perpendicular to the 
airplane center line) that contains the particular hound vortex being 
represented. Figure 6(a) shows a transverse section for the high-midwing 
configuration and figure 6(b) shows a plan view for the midwing 
conf iguration . 


This system of real and image vortices is only an approximation. 

It satisfies the condition of no flow across the fuselage boundary only 
at the transverse plane containing a bound vortex and its image and the 
transverse plane infinitely far behind the wing. In general, some flow 
will occur across the fuselage boundaries, and hence, to some degree at 
least, the corrective downwash velocities induced at the wing control 
points by the image vortices will be in error. 


Within these limitations it is therefore considered that wing- 
fuselage interference effects upon the wing spanwise airload distribu- 
tion can be obtained by: 


(1) Adding to the elements of the | 3 lJ matrix for the exposed wing 

a corrective matrix whose elements represent the vertical downwash 

velocities induced at the various wing horseshoe control points by the 
image vortex system within the fuselage boundaries. 


Each image vortex within the fuselage is of the same strength and 
sense of circulation as its real counterpart and there are as many image 
horseshoe vortices within the fuselage as are used to represent the wing. 


(2) Adding to the twist of the wing 


1 .e . 




a correction due 


to the vertical overvelocity field arising from the effects of fuselage 




is proportional to 


angle of attack. The correction applied to 

the quantity (o^ - a^^) where is the angle of attack of the wing 

root-section zero-lift line and a. is the incidence angle between the 


wing root-section zero-lift line and the fuselage center line . 


The above corrections are developed as follows. 
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The elements of the 
sidering the geometry of 


corrective matrix 
the following sketch; 


may be derived by con- 



The horseshoe vortex shown in this sketch is considered to be one 
of the image vortices, the location of which is given by equation (El) . 
The strength of circulation r of this image vortex i s eq ual to that 
of its real coimterpart . The plane of the horseshoe ooab«> is ^ an 
angle 7 to the horizontal plane Ogie and the bound vortex ab lies 
in the plane Omhg. The point k may be considered to be one of the 
downwash control points as described in appendix A. 

The total vertical downwash at the point k due to this image 

horseshoe vortex is the sum of the contributions due to each segment . By 
means of equation (A27) the contribution of the trailing vortex is 

- -r / i + cos ~ ^ ^ 

B yV 


B 


(E2) 
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that due to the trailing vortex hoo is 


boo li-K \ 


+ COS ^/ Sy' + H cos 7 


and that due to the bound vortex ab is 

(^os - cos 

D Ad 


r /c( 

w — = — (— 
ab kit \ 


(E5) 


(E4) 


where 


A = ck 
B = ^ 
D = fie 




(E5) 


From the geometry of the sketch the distances A, B, and D are found 
from 


A^ = ^Sy' + H cos 7^^ + ^s^' + H sin 7) 
B^ = ^Sy' - H cos 7)^ + ^82' + H sin 7^ 

2 2/ I \ 2 

D = s^' + /Sy' sin 7 - Sz cos 7J 

and the angles are obtained from 

D 




0' = tan ^3^= tan”^ 
fa 


(^' = tan”^-^^= tan”^ 
fb 


Sy' COS 7 + s^' sin 7 - H 


D 

s' cos 7 + s„' sin 7 + H 


A = tan ^ -^-= tan”^ 
ad 


|i = tan”^ tan ^ ^ 

- 


be 


(e6) 


/ (E7) 
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Summing the individual downwash' contributions (eqs. (E2) , (E3) , and 
(E4) ) gives the total downwash at k due to one image horseshoe 


''k’w 


+ H COS 7^ 1 + cos > 1 ^ ^ ^ cos - cos 6 ' )sjr 


A^ 


D 


^ Sy' - H cos 7^^1 + cos \ 




(E8) 


The downvash angles at the control points on the wing due to the 
fuselage images may then be written in matrix notation as 



(E9) 


in which the desired correction matrix for image vortex effects is js ^ , 

the elements of which are the locally applicable values of the function 
inside the brackets of equation (e 8) . The development for antisymmetrical 
conditions would be analogous to that for symmetrical conditions, the only 
change being that the sense of the vortices (and their images) on the 
left-hand side of the plane of symmetry in fig\ire 6 would be reversed. 


The second interference effect, that due to vertical overvelocity, 
may be taken into account as shown in the following analysis : With respect 

to the fuselage center line there is a transverse component of the free- 
stream velocity V^, given by (see sketch) 

= v(o^ - ai) 


Fuselage 



Root- section 
zero-lift line 


Since the fuselage displaces the air in this transverse flow field, there 
are local changes in this vertical flow velocity. 


If V 2 ' is equal to the total local vertical velocity in the pres- 
ence of the fuselage, then the overvelocity AV^ is 


AVz = Vz' - Vz 


(ElO) 
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This vertical overvelocity can he expressed as 


AV 


z " 



(Ell) 


Contour maps from which the value of AV2./V2 can be obtained as a 
function of vertical and lateral distances from the body center line are 
shown in figure 7 • 


These contoijrs were developed from the equations for uniform flow 
past an infinitely long circular cylinder from reference I7 as follows: 


V * P 

Z _ + a cos 29 

\ ~ r2 


(E12) 

(EI5) 


where a is the cylinder radius, R is the radial distance from the 
axis, and 0 is measured from a line normal to the axis and Vg* 
Thus 


Vz 


a^cos 20 


(E14) 


From equation (Ell) the angle-of-attack increment AXg at a control 
point due to the local overvelocity becomes 


Ax, 


g 




(EI5) 


Equation (EI5) in matrix notation is 



(E16) 


where is a diagonal matrix whose elements are the locally appli- 

cable values of AVg/V^; that is. 


— o 



(EI7) 
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pi] 

pa] 

pi] 

po] 

and the 


use of the matrices 

downwash matrix given by equation (A58) for the wing alone 
(the portion of the wing outboard of the fuselage) 

wing elasticity matrix given by equation (B5T) 
image-vortex matrix given by equation (E9) 

overvelocity matrix given by equation (EI 7 ) 
eqiaation 


fs} • [%]{i} 


the equation for wing equilibrium including fuselage -interference effects 
can be written as 


[Pi] * PiJ]« = |4qm<^|Jjat. + + agj- + [So]|aY - (EI 8 ) 


or, in a more convenient form, 






[I] . [4]1 { 1 } . [a,] 



(EI9) 


Equation (EI9) may be substituted without change for equation ( 21 ) for 
the wing load distribution. 
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APPENDIX F 

EQUATIONS FOR TAILLESS AND TAIL-BOOM AIRPLANE CONFIGURATIONS 


This appendix develops equations which, when properly inserted into 
the basic equations, allow solution for (l) a tailless airplane and 
(2) an airplane with the tail load entering the wing structure through 
tail booms. The equations are left in general form. Constant sweep 
angles and equal vortex spacing simplify the equations considerably. 


Tailless Airplane 

For the case of a tailless airplane longitudinal balance is us\ially 
accoiT5)lished by the deflection of controls on the wing. These control 
deflections alter the span load distribution over the wing while con- 
tributing to the balance of the airplane, and therefore terms expressed 
as functions of the control deflection 6 will appear in the lift- 
distribution equation (eq. (21 )) and the pitching -moment balance equa- 
tion (eq. (23)) in place of the P^ji terms. Since the lift produced by 

the control deflection is part of the wing lift, Pip in the lift balance 
equation (eq. (22)) will be zero. 

The development of expressions in terms of the control deflection 
follows . 

The required expression for 6 to be substituted into equation (21 ) 
consists of two ag coii5)onents; ag^, which is the apparent twist due to 

control deflection, and which is the twist due to section pitching 

moment with control deflected. From equation (C2) 



For class II twists due to section pitching-moment coefficient 
(type (f) in appendix C), the following equations may be written for 


{Mjcm and equation (Fl): 



(F2) 

(F3) 


88 


NACA TN 3050 


where 


h - 


1/2 

1 

1 

1 


0 

1/2 

1 

1 


0 

0 

1/2 

1 


0 

0 

0 

1/2 


For [m] and [t] in equation (F1), equations (B54) and (B55) 
written as 


where 


lo = 


W = - [jq] tin a] 

[t] = po] [cos a| 


(Fi)-) 

(F5) 


1/2 

1 

1 

1 


0 

1/2 

1 

1 


0 

0 

1/2 

1 


0 

0 

0 

1/2 


. 

• 

• 

• 



Substituting equations (F2), (F3)^ aJ^d (F5) into equation (F1) yields 
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or 


(S 




■ >i[5o][2ii] 


O O o Q o 

[^] + [^] [%] M{=-o} (f'S) 


Equation (f 6) gives the structural twist due to control deflection. If 
wind-timnel data are not available, theoretical expressions for 
in terms of the deflection & may be \ised. From reference l8^ 


H 


^ 1 

5 ■ 2 


sin 0^ - — sin 29r 




where 


0Q = sin 


-1 [ 2\||(1 - I ]) 

0Q = cos"l(2| - 1) 


I 

^ c 


and cf is the flap chord. Substituting these values for 0^ results in 

© 



” 0 

HiCVJ 

1 

II 

60 

T 

2\||(1 - 1 ) - 2\||(1 - |)(2| - 1) 


r— O — 1 


= -2 


'lid - 1)5 


{^) 


(FT) 


The ag^ term required can be obtained from the following expres- 
sion for the lift produced by control deflection (ref. 18); 


(°Z}6 = 2[^ - 0Q + sin 0^{6 


for a two-dimensional lift-curve slope of 2n. For a two-dimensional 
liftrcurve slope of iiIq^ this equation can be written as 




1 _ ^ + Qq 

" 7t 
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Substituting the values for 6 q results in 

o I— ° — 

{^2)5 = L“o] ^ cos-l(2| - 1) + i 2\1|(1 - 1) 

o 

Dividing both sides by [“o^l 


[5} (f 8 ) 


H' 


'HI 


1 _ i cos-l(2| - 1) + i 2(1(1 - 1) 


0 


(F 9 ) 


The total ag term to be substituted for the Pip term in the lift- 

distribution equation is obtained by combining eqviations (f 6), (F7), and 
(F 9 ), to give 


= 2q[lo] [2^ 

1 - i cos-1 (2| - 1) + i 2 ( 1(1 - I) 




(FIO) 



For the pitching-moment balance equation (25) the expression required 
to replace Pj may be obtained from the theoretical expression for flap 

pitching -moment coefficient (eq. (FT)) as 


PT = 2q 


2hc2 


o 


h(i - 1)5 (s) 


(Fll) 


Since in the system of equations given in matrix form (eqs. ( 2 l), 
(22), and (23)) there are two unknowns other than Pp, it wo-uld not be 
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possible to have more than one unknown in the expression to be substituted 
for Pijj in other words, the matrix in equations (FIO) and (F11) must- 

have a known distribution shape. Let 6^. be the unknown reference 
deflection} the matrix may then be written as where 

— > is any chosen distribution shape. For the reference deflection at 
the second station from the tip this equation becomes 


&l/Br 

1.0 

(^) = ^r ' S3/^r 


^root/^2 


For a constant deflection across the span, the distribution matrix 
becomes . 


Balancing Tail Load Entering Wing Through Tail Boom 

For the case where the tail load enters the wing structure through 
a tail boom, the distribution of load over the wing will be affected 
when the tail load changes and will vary in a different way depending on 

where the tail boom enters the wing. An expression to be substituted into 
the lift-distribution equation (2l) for developed in this sec- 

tion of the appendix. 

From figure 8, the beajnwise moment produced at the point P by the 
tail load entering the wing throiigh the tail boom can be written as 

Pm pm Pm ’ (Lj, 

Mp = — e^ sin + — yp cos Aj + -^ 2_ — 

Ten 2 2 ^ 2 cos An 


where Ap is the sweep angle of the elastic axis at the entry section 
of the tail boom. 
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In matrix form 


o 

1 


{m}pt = ^p2 ] M {It} + ^(®T sin A^ + yp cos At^[It} (F 12) 

where reference stations are defined as 


H 


0 

0 

1 
1 
1 
1 

Tip 


0 

0 

0 

1 

1 

1 


0 

0 

0 

0 

1 

1 


0 

0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

0 


>• — Tail-boom entry point 


Root 


[d] = 


0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

^e+1 


0 

0 

0 

0 

0 

^e+2 
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The equation for torque is 



- Aj - yp sin 


(FI 5 ) 


Substitution of equations (fIv), (F5), (F12), and (FI 3 ) into equa- 
tion (C2) gives an expression for the class II twist due to tail-boom 
entry: 



The column to be substituted into equation (21) is therefore 



This column will be multiplied by Pp in equation (2l). Equations (22) 
and ( 25 ) remain \onchanged for this configuration. 
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APPENDIX G 

METHOD OF REDUCING WIND-TUNNEL DATA 


Difficulty is usually experienced in applying data obtained from 
wind-tunnel tests to the design of a full-scale airplane, especially in 
those cases where of necessity the stiffness of the model differs from 
that of the full-scale airplane. The purpose of this appendix is to pre- 
sent a method of analysis by which model flexibility effects may be 
removed from the aerodynamic coefficients. 

Description of method .- The method utilizes equation (12), which is 

where 

H ' {W * W ^ W} 

When equation (12) is used for computing the lift distribution {ij 
for a given full-scale airplane of any given flexibility, the and 

rofl matrices used depend only on the particular configuration and the 

given flight conditions, except for certain aerodynamic -twist components 

of the matrix (see appendix C). There remain to be selected, then, 

applicable values for these components and for the JhIqJ matrix. 

These values may be determined with suitable data from wind-tunnel tests 
of a scaled model in conjvinction with equation (l2) as it applies to the 
model . 

The approach taken is usually applicable and is based on the assimip- 
tion that the following data are available from wind-tunnel tests of the 
model: U) the spanwise variation of the section normal-force coefficient 

obtained from integration of pressure data and (b) spanwise variation of 
section chord-line angle of attack with free-stream direction obtained from 
model deflection data. These data should be available for each of several 
root-section angles of attack and at each of several Mach n-umbers over the 
essentially linear range of section lift coefficient. 
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A downwash matrix for the model shovild be computed according 

to procedures of appendix A. In general the horseshoe system should be 
identical with that which will be used in determination of the airload 
distribution on the full-scale airplane. The matrix for the model [j^lJjyj 

will then be equal to the reciprocal of the model scale factor times the 
[Bl] matrix for the full-scale airplane, where the model scale factor is 
equal to the model span divided by the full-scale span. 


The final angle -of -attack matrix for the model may be written as 



■ {“) 


meas 



Ib 


(Gl) 


where 

{%•] 

Wib 


/a) 

V. Ji 


meas 


angle of attack of root section with undisturbed stream 
(measured) 

built-in twist (known) 


twist due to model flexibility (measinred) 


interference twist due to aerodynamic interference effects of 
neighboring bodies (fuselage, nacelles, external stores, etc.) 
.upon the wing (unknown) 

sum of measured or known values 


Equation (l2) for the model then becomes 



(G2) 


Since measured values of \}J are used, equation (G2) represents a 
system of independent equations (one for each reference spanwise station 
of the model) which do not require a simultaneous solution. 
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If assumed to be invariant with change in wing angle of 

attack, the following matrices may be computed from data taken at two 
different root-section angles of attack: 


H! - [“aleas - (“ijn 

M - {' J - {' j 


meas 


(G5) 

(G4) 


The matrices substituted for [^°'j]’nieas 

and in equation (G2) , since the term Wi. has been eliminated, 


and the applicable values of ^m^^ can be computed from 


= 


0 ^ 

1 


Uq Ax 




(G5) 


The interference twist 


may now be computed from equation (G2) 

by using the values of computed from eqmtion (G5) and the values of 

(ij and measured at any root-section angle of attack, for exam- 


, ^ meas 
pie |^Z2j and jo^j 


meas 


The foregoing procedvire indicates the simplest solution for the con- 
stants in the straight-line equations given by equations (G2) . If desired, 
equations (G2) may each be solved somewhat more accurately for K}. 

and ^ least-squares procedure which utilizes values of 

and /a} taken at several values of root-section angles of attack, 

u Jmeas 


If the interference twist a_ at any section is assumed to include 

Sib 

a component that varies with change in some section chord-line angle of 
attack, the method is still applicable to a close approximation. In this 

case the variable interference effects will appear implicitly in 

The approximation arises from the assumption (implied in the above 
solution) that for a body at some spanwise station j there exists an 
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interference twist at some other station i which is proportional to the 
change in geometric angle at i; whereas in actuality the interference 
twist at station i is proportional to change in geometric angle at the 
body station j . 

Justification for the assimiption can be shown as follows. Denote 
Interference twist by Oj, and the measured section chord-line angles at 

stations i and j by and . Interference twist at station i 

due to the body at station j is written as the sum of a constant plus 
a variable twist 

(■^)i ■ (%) * 

'i 

where k . . is the proportionality constant. The final angle of attack 
1 J 

at station i then is 

(Of). = ai + kijaj + 

and eqmtion (G2) for the spanwise station 

i [SlJ . [t] = + k. 

However, 

a. = a. + (a. 

J ^ \ J 

therefore^ 

i [SlJ ^ {j} - ”1 (^( 1 + ky ) '<■1 + ky (aj - ai) -f ( j (G«) 

The middle term is small compared to the first term and may be neglected 
since, for station i close to station j, the difference - a. is 
negligible, and for station i far from station j the interference 
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effect and therefore k is negligible. If the middle term is neglected, 
equation (g 6) becomes equivalent to equation (G2) and the factor 
1 + k,- is apparent in the values of mQ computed by equation (G5) 

even though it does not appear explicitly; that is, m^ computed by 
eqmtion (G5) is equal to the term m^^l + kj^j^ in equation (g6). 

Althou^ only interference twist was considered in the foregoing 
description, the method is obviously applicable for determining other 
aerodynamic twists (see appendix C), for example, that due to flap 
deflection. 

The foregoing method of analyzing data is relatively sinqple and 
straightforward and ,has the following advantages over other currently 
available methods: 

(1) It provides a means by which data obtained from wind-tunnel tests 
of a properly instrxunented elastic model wing can be evaluated for appli- 
cation to a full-scale wing of different elasticity distribution. 

(2) It evaluates the variation of effective section lift-c\irve slope 
and of aerodynamic twist across the span as influenced by the presence of 
the fuselage, the nacelles, and other bodies on or near the wing (these 
obviously include spoilers, ailerons, and flaps). 

(3) It determines the manner in which section lift-ciirve slope varies 
with Mach number. 

Although the variation of lift-curve slope with Mach mmiber can be 
obtained by the foregoing method by using wind-tunnel data at various 
Mach numbers, it is sometimes desired to determine the change (or changes) 
in lift distribution for section configurations for which suitable wind- 
tunnel data are not available, for example, variation of section lift 
with flap deflection c^^ . 

If the incompressible value of cj^ can be obtained or estimated 
by some means (e.g. , from tests for the flap deflected at some other sta- 
tion) and if the effective sweep angle at the new flap station is known, 
then compressible values of cj^ can be obtained by substituting the 

Cl values for the section lift-curve slopes in equation (A59) 

‘’5 


o r- 

\1 - M^cos^i^J4 

where uIq and m are the compressible and incompressible section lift- 
curve slopes, respectively. 
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It is therefore advantageous to evaluate the spanvise variation of 
The effective sweep angle for each section can be obtained from 

equation (A59) by using the previously meas\ired variation of iiiq with 

Mach number. For example, equation (A 59 ) can be rewritten as the linear 
equation (linear in squared terms) 

This equation can be solved for m and Am by a least-squares procedure 
if a series of equations are formed by substituting the values of m 
and M obtained for each of several Mach members. ° 


Alternatively, a solution for m and Am can be obtained by a 
graphical procedure (ref. 19 ). Plot the variation of mo with M in 
a fashio n such that t he abscissa (M-scale) is proportional to the fac- 

1//1 - M^cos^Am- The plotted points wi ll fall on a straight line 

which passes through the pole (m^ = 0, l/vi - m2cos2am = o), provided 
Am is correctly chosen and the law (eq. (A59)) applies. Note that 

t/ “ M^cos^Am ~ 1*0 when M = 0. A form is presented in figure 9 by 
which the mo values can be conveniently plotted on such an abscissa 
scale for each of a number of values of Am* For example, if a value 
of Am is tentatively chosen as Am = 50 °, draw a horizontal line inter- 
secting the right-hand ordinate at 30°. The intersections of this hori- 
zontal line with the Mach number lines are the abscissa locations for 
the indicated discrete values of M. Values of mo are then appro- 
priately plotted vertically above or below these intersections. The 
intersection of this horizontal line with the pole-distance curve gives 
the pole location in scale units to the left of the abscissa point M = 0. 
The desired value of Am is then the one which gives the most nearly 
linear variation of the plotted data with consideration given to the 
pole point . The desired m is the value of the intercept at M equal 
to zero. 


Illustrati ve example .- This section presents an analysis of wind- 
tunnel data obtained on a flexible model to determine the compressible 
section lift-curve slopes 

The model was a wing-f\iselage configuration with nacelles mounted 
below and forward of the wing on sweptforward struts attached to the wing 
semispan stations n = 0.57 and O. 65 . The wing had the following addi- 
tional pertinent physical characteristics; 


Aspect ratio 

Taper ratio 

Sweep at quarter-chord line, deg 


8.55 

o.4o 

55 


and the locus of aerodynamic centers was assumed to coincide with the 
quarter-chord line. 
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The model wing was instrumented with strain gages cemented to the 
surface of the wing steel spar along one semispan and with pressure 
orifices located along streamwise chord sections on the opposite semi- 
span at stations t] = 0.155^ 0 . 55 ^ 0*56, 0.75> aJid 0.92. 

Tests were made at root-section angles of attack of 0 ° and 6° for 
Mach numbers 0 . 50 , O.5O, O.7O, 0 . 75 ^ 8-nd 0.80. 


The spanwise variation in section angle of attack for each test con- 
dition was obtained from the root-section angles of attack in conjunction 
with strain-gage readings. The spanwise variation in section normal-force 
coefficient for each test condition was obtained by a spanwise fairing of 
the various local integrated chordwise pressure distributions. 

Elements of the matrix were computed for wing semispan stations 

11 = 0.10, 0.50, 0.50, 0.70, 0.85, 0.925^ and 0.975* resulting matrix 

is 



0.14869 
- 0.04084 
-0. 00287 
- 0. 00042 

-0.00015 
-0. 00008 
-0.00006 


0 . 04650 
0 . l 48 ll 
■0. 01260 
■0. 00068 
0.00017 
0. 00009 
■0.00006 


- 0. 01249 
-0. 05547 
0.07827 

-0 . 00426 
-0. 00055 
-0. 00025 
-0. 00015 


- 0. 00482 

-0. 00795 
-0.02713 
0.04582 
-0. 00522 
-0.00095 
-0. 00050 


-0. 00149 
-0. 0018a 
-0. 00285 
-0. 01109 
0 . o 4 i 53 
-0. 00617 
- 0. 00124 


-0. 00082 
-0.00094 
- 0. 00121 
-0. 00229 
-0. 01145 
0.05969 
-0. 00794 


-0. 00065 
-0. 00070 
-0 . 00084 
-0. 00150 
-0. 00287 
-0. 01507 
0.02917 


where the columns read down from the tip to the root and the rows read 
across from the tip to the root. For these same stations ACnc /4 and 
increments of section geometric angle of attack Ax due to the 6 ° change 
in root angle of attack were computed from data at each test Mach number 
and tabulated in table I. 


Substitution of 
equation (l 2 ) gives 


= ^ACnqcJ 


for 
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from which the values [®o]J w’ere computed and tabulated in table I. 

Values of itIq for t] = 0.975 are not shown since m^ is very sensitive 
to the manner in which the spanwise variation of ACj^ is faired in the 
region of the wing tip. These values for each wing station are plotted 
against Mach number in figure 10(a). 

Since the Pi] matrix is based upon the "wing alone" configuration, 

the body interference effects which vaiy with body angle of attack will 
be apparent in the values of lift-curve slope particularly for those 
stations near the wing root . 

It was assumed that the Prandtl-Glauert relationship was applicable 
for these data. Accordingly the values were plotted (by using a 

form similar to that shown in fig. 9) for various selected values of Aj^. 

The plots which yielded the most nearly linear variation of the data 
are reproduced in figure (l0(b) for each span station. In selecting the 
plot which gave the most linear variation of iiIq, consideration was given 

to the fact that at very low Mach numbers (approx. O. 5 O) the value of m^ 
could be affected by Reynolds number; whereas those obtained for high 
Mach numbers would reflect the effects of shocks, and so forth. These 
points were therefore given less weight in determining the best fit. 

Plots showing the spanwise variation of effective section sweep 
angle Aj^ and lift-curve slope m are given in figure 11. 
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APPENDIX H 

CALCULATION OF THE DIVERGENCE DYNAMIC PRESSURE 


AlthoTigh the sweptback wing is usually considered to be divergence - 
free, it is conceivable that, with a large external diverging torque 
such as may be contributed by a tip tank, the wing could diverge. This 
appendix is concerned with the calciilation of the dynamic pressure at 
which divergence will occur. 

An expression for the static lift on a flexible swept wing with an 
external tank has been derived in appendix D and, for a divergence inves- 
tigation, may be written as follows; 



o 


1 

kqmo 


f .] - H 



(HI) 


This is a matrix equation in which the elements in the column represent 
the lifts on the various spanwise segments of the wing. The total lift 
can be found from eqioation (HI) by multiplying the lift per unit span 

by the row matrix [2i^ where 2h defines the segment width in the 

spanwise direction. Then 


L 


= 2[_2hJ j^lj = 2[_2hl 




or 


:L2hJ 


- M 


-1 


H 


(H2) 


Equation (H2) represents the lift on a wing in static equilibrium. 
Under the condition of wing divergence the equilibrium wing lift would 
be infinite. 
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The divergence speed is then represented by the lowest value of q 
for which the lift in equation (H2) becomes infinite. For the right- 
hand side to be infinite, the determinant of the matrix whose inverse 
is given in equation (H2) must be zero, or 



S 




which can be written in the form 



0 


0 



1 

1 

rsii - sFai 

1 

rs-Ti - 

Sol 


4mo 

L y L j 

^“0 

Ly 

L^J 








(H3) 


The procedure is to solve for the lowest value of q which satisfies 
equation (H3). From the Cayley-Hamilton theorem, the dominant or highest 
modulus root A in the equation 



(H4) 


may be found by iterating the matrix In this case equation (H3) 

may be put in the form of equation (H4) by multiplying through by 


o 

1 ~ 

^nio 



The result is 
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Now 



therefore, 


iteration of the matrix product 


_ 0 


' -1 


0 




r 1 

Si 


S[A] 

1 

kmo 

Si 

+ 

S 2 ] 


(h6) 


will yield the lowest value of which in this case is the desired 
divergence dynamic pressure. 


The theory of matrix iteration may be found in reference 20; however, 
a simple outline of the steps required to iterate the matrix product (h 6) 
to obtain the divergence dynamic pressure is given here. First a trial 
column is chosen and this col\min is premiiltiplied by matrix (h 6) to 
obtain a result col\jmn. The elements of this resialt column are divided 
by the last element of the resvilt column and then become the elements of 
a second trial colimm (the last element will be imity, having been 
divided by itself ) . The second trial column is then premultiplied by 
matrix (h 6) to yield a second resiolt column. This procedixre is repeated 
until the same value is obtained for the last element in two successive 
res\ilt matrices. The reciprocal of this value is the desired divergence 
dynamic pressure q. 
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TABLE I.- MEASURED DATA AND RESULTING SECTION LIFT- CURVE SLOPES FOR 


WIND-TUNNEL MODEL OF APPENDIX G 



M = 0.50 

M = 0.50 

M = 0.70 

1 

Ac^cA 

Ax, radians 

irio 

Ac^cA 

Ax^ radians 

mo 

ACncA 

Ax, radians 

mo 

0.975 

.925 

.85 

.70 

.50 

.30 

.10 

15.897 
16 . 966 

20.775 

24.390 

30.490 

36.112 

37.485 

0.10135 

.10140 

.10148 

.10217 

.10317 

.10401 

.10455 

4.75 

5.41 

4.83 

5.84 

6.97 

7.19 

15.628 

17.460 

21.682 

27.023 

31.658 

56.112 

59.490 

0.09665 

.09672 

.09710 

.09849 

.10097 

.10294 

.10424 

5.08 

5.54 

5.98 

6.25 

6.69 

7.74 

14.558 

18.819 

22.998 

28.903 

55.669 

58.241 

45.587 

0.09117 

.09132 

.09192 

.09422 

.09852 

.10177 

.10592 

6.20 

6.08 

6.71 
6.79 
6.95 

8.72 


n 

M = 0.75 

M = 0.80 

ACncA 

Ax, radians 

mo 

AcncA 

Ax^ radians 

mo 

0.975 

.925 

.85 

.70 

.50 

.50 

.10 

15.510 

19.684 

23.814 

29.548 

34.857 

59.077 

45.157 

0.08972 

.08997 

.09066 

.09316 

.09772 

.10157 

.10382 

6.59 

6.46 

6.86 

7.15 

7.00 

9.09 

15.433 

20.631 

25.492 

51.804 

37.886 

42.650 

50.126 

0.08762 

.08792 

.08867 

.09152 

.09662 

.10082 

.10382 

6.97 

7.07 

7.45 

7.85 

7.59 

10.16 
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Figure 1 .- Angle-of -attack definitions and sign conventions. Positive 

angles shown. 
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Figure 2.- Typical vortex locations and location of pertinent points on 

an arbitrary wing plan form. 


110 


NACA TN 5050 



Figure 5.- Structural skeleton of outboard sections of left wing shown 

in plan view. 
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Figure 4 .- Plan view of left wing tip section. 
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Figure 5 *- Polar representation of 


left wing used in development of the 
matrix. 
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Airplane 


Locus for Rg= 


Fuselage boundary-^ 


Real vortex system 



'Wing plane 
at 3c/4 rearw/ard 
of leading edge 


-Corresponding 
image vortex system 


Point "k" at a distance of 
Sx' rearward of the transverse 
plane containing the bound vortex 
and its image and is on the locus ) 
of control points. 


(Note that a straight-line image bound 
vortex is used instead of its true shape, 
i.e., chord of the arc instead of the 
circular arc itself.) 


(a) Front view of typical vortex- image system (high-midwing configuration) . 
Figure 6 .- Diagrsim of vortex- image system for a wing-fuselage combination. 
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(b) Plan view of typical vortex- image system (midwing configuration) . 

Figure 6 .- Concluded. 


R/a = 2.00 +.01 +.02 +.03 +.04 
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R/a 




CTn 


Figure ?•- Concluded. 
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Figure 8.- Diagram of a representative wing — tail -boom combination showing 
entry of tail load on to wing structure (d is distance between stations 
measured perpendicular to streamline). 


Lift- curve slope, 


00 


M- scale for A =65° 



m 

Figure 9-- Plotting form based on equation mQ = u:- - 

Ml - M^cos^Aj^ 

graphical determination of effective section sweep angle A}^ 
a known variation of lift-curve slope with Mach number. 
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Spanwise 


C 

g 

o 


<D 

Q. 


o 

E 


O) 

Q. 

o 

a> 

O) 

> 

ZD 

0 

1 


c 

g 

t5 

CD 
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M- scale for A|^ = 35' 



5 .70 


0 

1 1 1 


0 

1 1 1 1 1 


0 .3 .5 .7 .75 .8 .85 

M- scale for A|^ = 35® 




.30 


0 .3 .5 .7 .75 .8 .85 

M” scale for A|y| = 58.5® 




tor A,^ =30^ 


(a) Linear scale. 


(b) Expanded M- scale. 


Figure 10.- Variation at several semispan stations of section lift-curve 
slope vith Mach number for linear scale and expanded M-scale. Wind- 
tunnel model of appendix G. 
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Spanwise station, t } 



Figure 11 .- Spanwise variation of effective sweep angle and section 

lift-c\irve slope m for wind-tunnel model of appendix G. 
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